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ABSTRACT 

The  two  principal  existing  methods  of  calculating  axially- 
symmetric  compressible  flow  in  turbomachines  are:    (1)  a  simplified 
one -dimensional  analysis,  and  (2)  numerical  methods  using  the  com- 
plete or  linearized  flow  equations.     The  first  is  not  satisfactory  for 
multi-stage  turbines  with  appreciable  wall  divergence;  the  second  is 
very  tedious  and  time  consuming.     The  purpose  of  this  investigation 
is  to  extend  the  approximate  methods,   successfully  used  in  calcu- 
lating incompressible  flow  in  compressors  with  constant  blade  height, 
to  the  analysis  of  compressible  flow  in  turbomachines  with  variable 
blade  height.    Assuming  that  the  blades  can  be  completely  defined  by 
the  exit  flow  angle,  and  neglecting  the  influence  of  downstream  blade3, 
the  analysis  is  made  considering  the  flow  between  successive  blade 
rows  only.     " IN  ith  these  restrictions,   subsonic  and  isentropic  super- 
sonic flow  patterns  can  be  determined  for  arbitrary  boundary  shapes 
as  long  as  separation  does  not  occur.    Average  losses  can  be  ac- 
counted for  by  the  use  of  a  polytropic  law,  and  the  effect  of  radial 
variations  in  stagnation  temperature  can  be  included  without  diffi- 
culty.     Examples  illustrating  the  flexibility  and  practical  value  of  the 
iteration  method,  and  the  rapid  convergence  of  successive  solutions 
are  given. 
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I.     INTRODUCTION 

In  the  modern  aircraft  gas  turbine,  the  flow  is  necessarily 
in  the  transonic  range  for  two  reasons:    (1)  to  keep  the  weight  and 
size  of  the  machine  at  a  minimum,  and  (2)  to  fit  the  mass  flow  and 
pressure  ratio  of  the  compressor.     Earlier  turbojet  engines  had  single 
stage  transonic  turbines,  but  present  engines  employ  multistage  tur- 
bines with  rapid  channel  divergence  in  order  to  use  effectively  the 
increased  pressure  ratios  and  keep  the  axial  velocities  compatible 
with  the  blade  speed.     In  these  types  of  machines,  the  flow  at  the 
mean  blade  height  is  mainly  asdal,  and  the  aspect  ratio  of  the  blades 
is  usually  above  unity  for  weight  saving  considerations.    As  long  as 
the  channel  wall  divergence  is  not  too  great  so  that  separation  can 
be  avoided,  there  is  no  reason  for  low  efficiency,  particularly  in 
view  of  the  favorable  pressure  gradient  through  the  turbine. 

Two  distinct  design  problems  occur  in  practice  which  Marble 
(Reference  15)  has  classified  as  the  direct  and  inverse  problems  of 
turbomachine  design  in  analogy  with  the  corresponding  classical 
problems  in  the  theory  of  finite  wings.   The  direct,  or  "off-design" 
problem  occurs  when  the  blade  speed,  blade  shape,  boundary  configu- 
ration, and  fluid  state  ahead  of  all  blades  is  given,  and  the  three- 
dimensional  velocity  field,  blade  loading,  and  distribution  of  energy 
are  to  be  determined.     The  inverse  problem  occurs  in  the  original 
design  of  turbomachine ry,  when  the  blade  loading,  blade  speed, 
boundary  configuration,  and  fluid  state  ahead  of  all  blades  is  pre- 
scribed, and  the  three-dimensional  velocity  field,  blade  shape,  and 
distribution  of  energy  in  the  field  are  to  be  determined. 

The  exact  analysis  of  viscous,  compressible  flow  through  an 
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axial  turbomachinc  will  probably  never  be  possible  because  of  the 
non-linear  character  of  the  equations  governing  the  flow.     However, 
Via  design  engineer  has  a  great  need  for  an  analysis  suitable  for 
engineering  purposes,  and  may  have  no  u3e  for  a  lengthy,  coxr-pli- 
cated  and  costly  procedure  even  if  it  were  available.     He  needs  an 
extremely  flexible  method  of  analysis  which,  while  being  as  simpli- 
fied as  practicable,  still  produces  information  accurate  enough  for 
effective  design  work. 

The  equations  governing  three-dimensional,  compressible 
channel  flow  are  greatly  simplified  if  terms  involving  viscosity  are 
neglected.     In  view  of  the  fact  that  viscous  effects  actually  are  very 
slight  except  near  the  boundaries,  and  as  long  a3  separation  is  avoided, 
the  assumption  of  zero  viscosity  seems  justified.     Viscosity  effects 
could  be  accounted  for  as  an  added  correction  to  the  basic  solution. 
However,  with  the  present  day  turbines  which  operate  in  the  high- 
subsonic  range,  it  is  obvious  that  compressibility  effects  can  not 
be  neglected. 

One  of  the  principal  methods  available  to  the  design  engineer 
for  the  analysis  of  axially- symmetric,  compressible  flow  in  turbo- 
machines  is  a  simplified  one -dimensional  analysis.     This  method  is 
used  extensively  in  preliminary  design  considerations  and  is  simple 
and  straightforward.   Since  only  the  exit  flow  angles  from  the  blade 
rows  at  the  mean  blade  height  are  considered,  the  radial  distributions 
of  velocity  and  energy  are  not  determined.    Although  this  analysis 
gives  reasonably  accurate  results  for  a  single  stage  turbine  with  a 
fairly  high  hub  ratio  and  small  channel  divergence,  the  one -dimensional 
method  is  inadequate  for  large  channel  divergence  where  the  radial 
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distribution  of  velocities  becomes  significant.     It  is  particularly  in- 
adequate for  the  design  of  multistage  turbines  with  channel  wall  di- 
vergence since  it  is  difficult  to  get  the  proper  balance  of  work  on  the 
individual  stages.     This  work  balance  is  strongly  affected  by  the  radial 
distribution  of  the  axial  and  radial  velocities.     Since  the  one-dimen- 
sional analysis  is  done  at  the  mean  blade  height,  it  gives  no  informa- 
tion as  to  the  radial  distribution  of  velocities. 

The  second  principal  method  of  analysis  of  axially- symmetric 
compressible  flow  in  turbomachines  is  by  numerical  methods  using  the 
complete  or  linearized  flow  equations  to  solve  for  the  complete  velocity 
and  energy  field  through  the  machine.     Since  the  determination  of  the 
entire  flow  field  involves  six  unknowns  at  each  point  in  the  flow  (the 
three  components  of  velocity,  the  pressure,  the  temperature,  and  the 
density),  the  analysis  using  these  numerical  methods  is  necessarily 
lengthy  and  extremely  complicated  and  cannot  be  considered  suitable 
from  the  point  of  view  of  the  design  engineer.     Monroe  (Reference  7) 
formulated  the  idealized  problem  in  terms  of  a  stream  function  for  the 
velocities  in  the  meridional  plane,  and  solved  the  resulting  non-linear 
differential  equation  by  a  simultaneous  application  of  the  relaxation 
technique  of  Southwell  and  an  iteration  process.     Wu  (Reference  9)  has 
utilized  various  relaxation  methods,   matrix  systems,  and  finite  differ- 
ence schemes  to  solve  the  complete  hydrodynamical  equations  for  the 
entire  flow  field  in  many  types  of  problems.     Vavra  (Reference  8), 
Reissner  (Reference  10),   Goldstein  (Reference  11),   Vvislicenus  (Ref- 
erence 13),  and  many  others  have  employed  various  techniques  to  solve 
the  complete,  or  nearly  complete,  equations  by  some  numerical  method. 
Although  these  numerical  methods  are  not  suitable  for  design  purposes, 
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they  are  quite  general  and  can  be  used  for  centrifugal,  mixed  flow  or 
axial  flow  machines  with  blades  of  any  aspect  ratio. 

In  an  effort  to  simplify  the  problem  of  solving  for  the  complete 
flow  field,  a  linearized  analysis  has  been  provided  by  Marble  (Refer- 
ences 5  and  15)  and  Fabri  (Reference  6).     Marble  (Reference  15)  pro- 
posed a  linearized  analysis  of  the  problem,  satisfactory  for  high  aspect 
ratio  blades,  to  allow  an  approximate  treatment  of  the  general  blade  row 
with  prescribed  wall  geometry.     He  assumed  that  the  trailing  vorticity 
could  be  considered  to  be  transported  downstream  by  the  mean  axial 
flow  and  not  by  the  perturbation  velocities,  that  is,  disturbances  in  the 
radial  and  axial  velocities  are  small  in  comparison  with  the  mean  axial 
velocity,  and  that  the  blade  row  is  made  up  of  an  infinite  number  of 
blades  of  finite  chord.   The  analysis  for  the  inverse  problem  was  carried 
out  in  detail  and  was  found  to  allow  a  general  and  reasonably  simple  so- 
lution.    In  Reference  5,  Marble  extended  the  analysis  to  the  direct  prob- 
lem of  off-design  operation,  and  to  the  study  of  mutual  interference  of 
neighboring  blade  rows  in  a  multistage  axial  turbomachine.     He  provided 
examples  of  axial  and  conical  flow,  and  stated  that  the  simple  linearized 
solution  was  sufficiently  accurate  if  the  vorticity  effects  were  not  large. 
A  second  order  linearization  was  given  to  handle  problems  with  greater 
vorticity  effects.     However,  the  computation  required  to  find  the  three- 
dimensional  flow  in  any  particular  case  is  rather  lengthy  using  the  lin- 
earized solutions,  and  involves  some  type  of  numerical  integration.    A 
simple  exponential  approximation  was  then  introduced  to  simplify  the 
computation,  and  the  results  using  this  approximation  compared  favor- 
ably with  the  detailed  results  of  the  linearized  solution. 

The  analysis  of  flow  in  turbomachine ry  where  the  channel  bound- 
aries are  purely  cylindrical,  or  vary  only  slightly  from  a  cylindrical 
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shape,  has  been  extensively  investigated  by  Rannie  (References  1  and  14), 
Traupel  (References  2  and  3),  Sinnette  (References  4  and  16),  and  others. 
In  these  analyses,  the  radial  momentum  equation  used  in  the  treatment  is 
of  the  form, 

e    dr  r  (1), 

where  the  radial  force  and  other  momentum  terms  have  been  neglected. 
The  distribution  of  axial  velocity  is  approximated  by  assuming  axial  sym- 
metry, and  since  the  centrifugal  force  within  the  rotating  fluid  body  is  bal- 
anced only  by  the  radial  pressure  gradient,  the  resulting  flow  calculated 
in  this  manner  is  in  reality  that  which  must  exist  far  downstream  of  the 
blades  where  radial  velocities  and  accelerations  have  vanished.     No  in- 
formation is  provided  by  these  analyses  on  how  rapidly  the  change  in  ve- 
locity pattern  takes  place  as  the  fluid  passes  through  the  blade  row. 

The  purpose  of  tide  thesis  is  to  present  an  approximate  method 
for  the  analysis  of  compressible  flow  in  modern  axial  turbomachines  which 
gives  sufficiently  accurate  information  for  design  purposes  and  yet  retains 
the  simplicity  necessary  to  be  of  practical  value  to  the  design  engineer. 
This  is  the  logical  step  forward  from  the  one -dimensional  method  to  the 
analysis  of  compressible  flow  in  axial  turbomachines.     Methods  which 
have  been  successfully  employed  in  the  incompressible  case  (Reference  1), 
have  been  modified  and  extended  to  the  analysis  of  compressible  flow  in 
arbitrarily  shaped  axially- symmetric  channels.     It  is  obvious  that  little 
simplification  over  the  so-called  "numerical"  methods  can  be  obtained  if 
the  complete  flow  and  energy  fields  are  to  be  determined  in  the  analysis. 
In  this  thesis,  the  flow  is  analyzed  at  the  station  immediately  downstream 
of  the  blade  row  by  a  consideration  of  the  upstream  conditions,  channel  con- 
figuration, and  blade  shapes.     The  blade  trailing  edge  angle  only  is  consid- 
ered, and  the  complete  flow  field  through  the  blade  row  is  not  determined. 
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The  method  of  analysis  presented  is  extremely  flexible  and  has 
been  used  to  solve  the  direct  and  inverse  problems,   subsonic  or  super- 
sonic, for  arbitrary  axial  boundary  configurations.     In  the  solution  of 
the  direct  problem,  the  choice  of  blade  speed,  blade  spouting  angle, 
boundary  configuration,  and  upstream  conditions  can  be  arbitrary,  and 
the  three  -dimensional  velocity  field,  blade  loading,  and  distribution  of 
energy  in  the  field  are  determined.    In  the  inverse  problem,  for  a  given 
blade  loading,  blade  speed,  boundary  configuration,  and  fluid  state  ahead 
of  all  blades,  the  blade  spouting  angle,  velocity  field  and  energy  distri- 
bution can  be  determined.     By  a  simple  change  of  3ign  in  the  equations 
used,  the  method  can  be  applied  equally  well  to  an  axial  compressor  as 
to  an  axial  turbine.     Although  this  method  is  not  as  general  as  those  U3ed 
by  Wu  (Reference  9),  for  example,  it  is  more  flexible  and  practical  for 
blading  of  aspect  ratio  greater  than  unity,  and  is  sufficiently  accurate 
in  view  of  the  assumptions  used  in  the  treatment. 

The  analysis  is  carried  out  on  the  basis  of  steady,  axially- 
symmetric  flow  with  infinitely  many  blades.     The  fluid  is  assumed  to 
be  compressible,  but  all  other  real  fluid  effects  are  neglected  except 
as  they  may  be  approximated  by  a  "polytropic  efficiency1 '  and  the  use 
of  the  polytropic  exponent  "n".    As  a  result  of  this  assumption,  the  pos- 
sibility of  separation  of  flow  at  the  channel  walls  or  along  the  blades  is 
excluded.     Boundary  layer  effects  are  also  not  considered  in  the  anal- 
ysis.    In  the  supersonic  solution,  the  occurrence  of  shock  waves  is 
excluded.     Radial  or  nearly  radial  blades  are  presupposed  so  that  the 
radial  force  can  be  neglected,  and  the  radial  momentum  equation  used 
in  the  treatment  has  the  form, 

>iR     .   }d_a9u        da  (2) 
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This  assumption  seems  justified  in  view  of  the  investigations  of  -\arlsson 
(Reference  19).  who,  by  a  two-dimensional  analysis  of  incompressible 
flow  with  an  infinite  number  of  blades,  showed  that  the  effect  of  the  radial 
force  resulting  from  normal  blade  twist  is  of  negligible  magnitude. 

The  effect  of  downstream  blades  on  conditions  upstream  is  ne- 
glected, and  the  blade  shape  is  specified  by  the  trailing  edge  angle  only. 
Simply  stated,  a  blade  row  is  assumed  whose  only  effect  on  the  flow  is 
to  turn  it  through  a  specified  angle  distribution  which  is  a  function  of 
radius  across  the  channel.    An  assumption  similar  to  this  has  been  used 
successfully  in  the  analysis  of  flow  through  a  compressor  (References 
1  and  14),  and  should  be  even  more  applicable  to  turbines  where  the  blade 
solidity  is  high*     The  equations  are  particularly  suited  to  the  analysis 
of  flows  with  a  prescribed  radial  total  enthalpy  distribution. 

Subject  to  the  above  assumptions,  the  resulting  u exact"  equations 
are  solved  by  an  approximate  method  using  a  simple  iteration  procedure. 
An  integral  for  the  axial  velocity  is  derived  in  terms  of  Known  upstream 
conditions  and  the  downstream  pressure.   The  pressure  distribution  is 
found  by  integrating  the  radial  momentum  equation  and  applying  the 
energy  equation  at  a  boundary.     The  maas  flow  equation  is  used  to  solve 
for  the  constant  obtained  from  the  axial  velocity  integral.     The  iteration 
is  started  by  assuming  an  axial  velocity  distribution  and  continued  until 
the  solution  is  reached,  i.  e. ,  until  the  velocity  distribution  produced  by 
an  iteration  is  the  same  as  that  assumed  for  that  particular  iteration. 
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II.     FUNDAMENTAL  EQUATIONS 

A.     Development  of  Basic  Differential  Equations 

The  differential  equations  are  developed  assuming  steady, 
compressible  flow  of  a  fluid  through  an  axially  symmetric  cylindrical 
channel  of  arbitrary  shape.    An  infinite  number  of  blades  is  assumed 
so  that  the  velocity  components  are  considered  independent  of  angle 
about  the  axis  of  the  machine.     The  fluid  is  assumed  perfect  with  no 
losses  of  any  sort  except  as  these  losses  can  be  approximated  by  a 
"polytropic  efficiency"  and  the  use  of  the  polytropic  exponent  "n". 
All  boundary  layer  and  viscosity  effects  are  also  neglected. 

In  order  to  generalise  the  equations  for  an  arbitrary  (but 
reasonable)  channel  shape,  it  is  necessary  to  consider  a  special  type 
of  blading.     These  blades,  infinite  in  number,  have  no  other  effect  on 
the  flow  than  to  turn  it  through  a  prescribed  angle  distribution.     Thus 
the  effect  of  downstream  blades  on  upstream  conditions  is  neglected. 
Attention  is  directed  to  the  flow  patterns  immediately  before  and  after 
the  blade  rows  for  either  the  stator  or  the  rotor. 

In  Figure  i,   a  stream  annulus  of  small  radial  extent  is  shown, 
and  absolute  velocity  components  are  considered  at  stations  1,  2,  and 
3.     Generally,  there  will  be  a  shift  of  the  stream  surface  in  passing 
through  either  stator  or  rotor  blading,  and  three  components  of 
velocity  will  be  present  at  each  station. 

With  the  above  assumptions,  the  continuity  equation  is  written 
in  differential  form  as: 


e, %  r, <*-,  '  fzvfzrzd^  =  e3 *3 r, dr3 
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FIG.  I—  CHANNEL     AND     STREAM     ANNULUS 

For  a  row  of  stator  blades,  the  energy  equation  must  hold  along 
a  stream  surface,  and  for  a  row  of  rotor  blades,  the  energy  equation 
must  hold  along  a  relative  stream  surface.     Hence,  in  terms  of 
absolute  velocities  at  stations  1,  2,  and  3,  we  have: 
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Neglecting  the  radial  force,  the  radial  momentum  equation 
hold  at  all  stations.    In  absolute  velocities  and  omitting  subscripts, 
:  3  equation  is: 


_L  22L    = 
f  dr 


(2). 


If  the  flow  is  considered  to  be  isentropic  along  a  relative  streatr 
surface,  the  following  relations  must  hold  along  this  surface: 
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If  a  polytropic  flow  is  considered,     ^      will  be  replaced  by    n    in 
Equations  6  and  7.     The  above  equations  are  fundamental  in  the  solution 
of  the  direct  an  I  inverse  problems . 

B.     Transformation  to  Non-dimensional  Foriri 


For  ease  in  handling  the  equations,  and  to  avoid  the  use  of  large 
numoers  in  calculation,  the  derivation  is  completed  in  non-dimensional 
form.     The  notation  used  is  given  in  Appendix  A.     Axial  and  radial  co- 
ordinates are  made  non-dimensional  by  referring  them  to  a  standard 
dimension,  usually  the  outer  radius  of  the  annulus  at  station  i.    Velo- 
cities are  made  dimensionless  by  referring  them  to  a  standard  tip  speed. 
The  following  notation  is  used: 

It  -  outer  radius  of  the  annulus  at  station  1. 

o 

cj  *\  ^  -  reference  blade  speed. 

o  =   ^  ft  -  dimensionless  radial  component  of  velocity. 

A=       „  -  dimensionless  tangential  component  of  velocity. 

-  dimensionless  axial  component  of  velocity. 

-  dimensionless  pressure. 

-  dimensionless  density. 

-  dimensionless  temperature. 

-  gas  equation. 
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777   =    — ^-^ —  -  dimensionless  mass  flow, 

/   =   — —  -  dimensionless  radial  distance. 

7  Ro 

£  =  -Z—  -  dimensionless  axial  distance. 

With  this  notation,  Equations  (3),  (4),  (5),  (2),  (6),  and  (7)  become 
respectively: 

<?,  *,  f,  *i      -   <?e  %  f.  «%  =    tfs  *,  f,*t,  (8> 

£**-■£$  + 1«+  *'+  f)  -  %  i+  r  (<+  *>*>  £rH  (9) 

-klf  -  i -•#-*$  (n) 

If  -  ft/-  (IT  (12) 


<?*   _   //>*)*_  /_?j 


<?£  l^/  ** 


L  ' 

r 
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C .    Equations  for  the  Direct  Problem  (or  an  Arbitrary  Channel  Shape 

The  basic  integral  for  the  axial  velocity  will  be  derived  for  the 
rotor  (stations  2  to  3)  for. the  isentropic  case  and  simplified  to  the 
integral  for  the  stator.     The  complete  derivation  is  given  in  Appendix 
B  for  the  isentropic  direct  problem  and  in  Appendix  C  for  the  polytropic 
direct  problem. 

Differentiating  the  left  side  of  Equation  (10)  in  respect  to  the 
radius         /         and  the  right  side  in  respect  to      £        gives: 
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JV  dtkt_  ctteJehs  m   \JL  *  (3)  ,  >  <Llg\  A \ <£)+  4rA\ctf     ( 1 4) 

Using  Equations  (8)  and  (13),  differentiation  of  the  first  term  of  the  right 

side  of  the  above  equations  results  in: 

j_ 
JLd-fJk)  -    J-  *3  +  JL  £*  A  /Ar) hi3  { 1 5) 

Substituting  Equations  (11)  and  (15)  into  (14),  and  using  the  differential 
form  of  the  mass  flow  relation  (i.Cquation  (8)  )  for  the  ratio  g  , 

the  following  differential  equation  is  obtained: 

2  -*■  f 


(16). 
d&3 


For  a  cascade  with  an  infinite  number 'of  blades,  the  exit  flow 
is  parallel  to  the  tangent  of  the  camber  line  at  the  trailing  edge,  and 
its  direction  is  independent  of  the  direction  of  flow  approaching  the 
cascade  .    For  most  practical  applications ,  the  direction  of  flow 
leaving  the  cascade  is  very  nearly  constant  and  independent  of  the 
inlet  direction  throughout  the  design  operating  range  .     Vi  ith  this 
assumption  for  the  rotor,  it  is  therefore  assumed  that  the  rotor 
spouting  angle       fl3        is  a  known  function  of  radius  at  station  3.    Using 
the  sign  convention  shown  in  Figure  8,  it  is  therefore  possible  to 
describe  the  relation  between         r\3        and      <p3      along  the  radius  at 
station  3  with  a  function  G  ( f3)  such  that: 
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Eiirninating      ^3    from  the  differential  Equation  (16)  by  use  of 
Equation  (17)  and  integrating,  the  following  integral  for  the  rotor  axial 
velocity  for  the  direct  problem  is  obtained: 


4 


+[*MJ  ft  fat, 
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It  is  to  be  noted  that  no  linearisation  has  been  performed  in 
arriving  at  the  above  integral.     The  expression  for    $~    is  quite  compli- 
cated,  but  its  numerical  evaluation  by  a  simple  iteration  process  to  be 
described  later  is  straightforward.     Once  a  set  of  streamlines  have 
been  assumed,   all  factors  involving  the  upstream  conditions  become 
known  functions  of     f3       .     In  the  solution  of  the  direct  problem,  it 
is  assumed  that    G    (f3)     is  a  given  function  of    £3     .     Thus  in  the 
integrand,  P,  and        ,  > —       are  the  unknowns  which  must  be  deter  mined 
in  the  iteration  process  in  order  to  evaluate  the  integral. 

An  alternate  form  for  the  square  bracket  in  the  integrand  can  be 
obtained  by  using  the  energy  and  radial  momentum  equations  at  station  2: 

(19). 


In  the  case  that  P^*  =  1  (isentropic  expansion  to  station  2  from  stagna- 
tion conditions),  the  bracketed  term  in  the  integrand  may  be  further 
simplified  to: 

(20). 


fcf  ('-&-&'•*$ 
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The  constant  of  integration  C,  is  determined  by  integrating 
Equation  (8).  The  integrated  form  of  this  relation  equates  the  total 
flow  rate  upstream  to  that  at  station  3: 

771     -f^MY^s^  (2D. 

Using  Equation  (13),  the  integral  form  of  the  radial  momentum 


equation  becomes: 


p3~ 
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The  constant  of  integration  C,  is  obtained  by  applying  the  energy 
relation  (Equation  (10)  )  at  a  boundary. 

Having  assumed  a  set  of  streamsurfaces  for  a  particular 
boundary  configuration,  the  value  of       ,  ,  can  be  approximated 

from  the  values  of  the  slope  and  rate  of  change  of  slope  of  these  stream- 
surfaces  at  station  3.     At  each  point  on  a  streamline  at  station  3,  the 
following  relation  must  hold: 

where        (    ,  ^  J        is  the  slope  of  the  streamline  in  question.     Thus 
the  value  of         ,  „*       becomes: 

In  order  to  approximate  the  value  of        ,  J         at  points  across  the 
channel,  it  is  therefore  necessary  to  estimate  the  slope  and  curvature 
of  the  assumed  streamlines  at  each  point. 
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For  the  solution  of  the  direct  problem  for  the  rotor  with  an 
arbitrary  channel  shape,  Equations  (10),  (13),  (17),  (18),  (21),  (22), 
and  (24)   are  applicable. 

The  derivation  of  the  integral  for  the  axial  velocity  in  the  case 
of  the  stator  is  carried  out  in  a  similar  fashion  using  the  form  of  the 
energy  relation  given  as  Equation  (9).     Using  the  sign  convention  of 
Figure  8,  Equation  (17)  for  the  stator  (stations  1  to  2)  becomes: 

*a  =  F(fz)  <Pa  (25). 

The  integral  for    9-   obviously  takes  the  same  form  as 
Equation  (13),  and  is  given  here  only  for  completeness: 

_  Jk  &* 


1     f£i(  &  \jL  &* 
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jr  the  solution  of  the  direct  problem  for  the  stator  with  an 
arbitrary  channel  shape,  Equations  (9),  (12)„  (21),  (22),  (24),  (25), 
and  (26)  (with  appropriate  subscripts)  are  applicable. 

D.     Equations  for  the  Inverse  Problem  for  an  Arbitrary  Channel  Shape 

In  the  inverse  problem  for  the  rotor,  the  same  assumptions  and 
basic  differential  equations  are  used  as  in  the  direct  problem  with  the 
exception  that  now  the  blade  relation,  Equation  (17),  is  replaced  by 
the  following  given  function  which  describes  the  work  output  of  the 
rotor  as  a  function  of  the  radius  at  station  3: 


"(f,)°  AV£a, 


(27). 
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The  function  K  (f  )  gives  the  dimensionless  rate  of  change 

of  angular  momentum  per  unit  mass  flow  at  the  particular  radius,  and 
represents  the  energy  taken  out  of  the  stream.     By  multiplying  this 
function  by  suitable  constants  it  can  readily  be  converted  to  B.  T.  U. 
per  pound  of  mass  flow. 

The  complete  derivation  of  the  axial  velocity  integral  is  given 
in  Appendix  D  for  the  isentropic  case  and  in  Appendix  E  for  the  poly- 
tropic  process.     For  the  isentropic  inverse  problem,  the  derivation 
remains  the  same  up  to  and  including  Equation  (16),  -which  may  be 
■written  in  a  slightly  modified  form  as: 


fi 


+  <t>d*s 


*Jf.  (28). 


Using  Equation  (27)  to  eliminate         A-  ,  the  integral  is 


3 


found  to  be: 

at  4 


*3  \         %z  M& 


dft  +C, 

(29). 


Again  the  numerical  evaluation  of  this  integral  by  the  iteration 

process  to  be  described  later  is  straightforward  once  a  set  of  stream- 
ed O3 

lines  has  been  assumed.     The  unknown    P,  and  M         appear  in 

cL  f3 

the  integrand  as  in  the  direct  problem.     Because  of  the  nature  of  t 
differential  equation,  however,  the  axial  velocity    < :■-    aiao  appears  in 
the  integrand. 
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Since  the  mass  flow  and  radial  momentum  integrals  remain 
unchanged,  Equations  (10),  (13),  (21),  (22),  (24),  (27),  and  (29)  are 
applicable  to  the  solution  of  the  isentropic  inverse  problem. 

If  the  inverse  problem  could  be  conceivably  applied  to  a  stator 
it  would  consist  of  prescribing  values  of       A2     as  a  function  of    <f2      , 
and  solving  for  the  resulting  velocity  field,  distribution  of  energy  and 
blade  shape  (as  described  by  the  blade  trailing  edge  angle).    In  this 
case,  the  integral  for  the  axial  velocity  simplifies  to: 


d£  +  c, 


Jfi  (30) 

where      ^ztfz)       *s  a  prescribed  function  of  the  radius  at  station  2. 
Equations  (9),  (12),   (21),  (22),   (24),   and  (30)  (with  appropriate 
subscripts)  would  be  applicable. 

E .     Changes  in  Equations  for  the  Compressor  Solution 

In  the  development  of  the  equations  in  the  preceding  sections, 
the  turbine  sign  convention  given  in  Figure  8  has  been  adopted.    If 
the  equations  are  to  be  used  for  compressor  design,  Equation  (17)  for 
the  direct  rotor  problem  becomes,  in  accordance  with  Figure  9: 

*j  =  f.  ~  &(&  +>  (31)- 

and  the  energy  relation,  Equation  (10),  is  changed  to: 
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Using  these  equations ,  the  derivation  is  carried  out  in  exactly  the 
same  manner  as  set  forth  in  the  preceding  sections  for  the  turbine. 

For  the  inverse  problem  for  the  compressor ,  the  work  out- 
put function  (Equation  (27)  )  becomes: 

"*>-**•-<*  (33) 

Using  Equations  (32)  and  (33)  in  place  of  Equations  (10)  and  (27),  the 
relations  for  the  inverse  compressor  problem  are  derived  in  the  same 
manner  as  for  the  turbine. 
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III.     PROPOSED  ITERATION  PROCEDURE 
A.     Procedure  for  the  Direct  Problem 

The  general  scheme  ol  iteration  is  exactly  the  same  for  the 
direct  and  inverse  problems,   subsonic  or  supersonic,  isentropic  or 
polytropic.    In  Appendix  F,  one  complete  iteration  is  given  for  an 
isentropic  supersonic  solution  (Example  Vb)  and  a  polytropic  subsonic 
solution  (Example  VIII  b)  for  the  direct  rotor  problem.     The  stator  solu- 
tion is  obtained  by  the  same  procedure  using  the  appropriate  axial 
velocity  integral  and  energy  equation. 

In  order  to  understand  the  technique  of  iteration  for  the  direct 
problem  described  below,  the  reader  is  urged  to  refer  frequently  to 
Examples  Vb  and  VIII b  of  Appendix  F.     For  an  arbitrary  channel 
shape,   seven  equations  are  used  in  the  solution.     For  the  rotor,  these 
relations  are  Equations  (10),   (13),  (17),  (18),   (21),   (22),  and  (24).     In 
the  examples  worked  out  in  Appendix  F,  the  functions  "/  \$z}  a**d  ^  KW 
have  been  assumed  to  be  constant  and  independent  of  radius  in  order 
to  simplify  the  axial  velocity  integral  and  the  computation.     This 
assumption  is  by  no  ineans  necessary.     For  arbitrary  given  functions, 
F  ( £  J         and  G  (43)     »  the  exponential  integrals  of  Equation  (18) 
should  be  evaluated  at  various  stations  across  the  channel  by  plotting 
the  integrand  and  determining  the  value  of  the  integral  graphically. 

In  the  solution  of  the  direct  problem,  it  is  assumed  that  the 
boundary  configuration,   blade  speed,  blade  trailing  edge  angle,  and 
upstream  conditions  are  prescribed.     It  is  further  assumed  that  the 
mass  flow  rate  is  fixed  and  is  compatible  with  the  given  upstream 
conditions. 
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Once  the  problem  has  been  set  up,  the  first  step  is  to  assume 
a  set  of  streamsurfaces  as  shown  in  Figures  3,   4,   and  5.     For  most 
applications,  these  streamsurfaces  should  be  chosen  on  the  basis  of 
equal  annular  areas  in  a  plane  perpendicular  to  the  axis  of  the  machine. 
The  accuracy  of  the  solution  can  be  improved  by  increasing  the  number 
of  stations  across  the  channel.     However,  it  is  believed  that  for  most 
purposes,    six  or  seven  steps  will  give  sufficient  accuracy  for  design 
work. 

Having  assumed  a  reasonable  set  of  streamlines,  the  second 
step  is  to  estimate  the  axial  velocity  distribution  at  the  downstream 
station.     Taking  average  or  middle -channel  values  for  the  prescribed 
upstream  conditions,  the  mass  flow  relation  and  the  energy  equation 
can  be  combined  into  one  equation  in  terms  of  <j>,.     Using  average 
values  in  Equations  (10)  and  (21),  the  relation  is  obtained  as  follows: 


r  -  Z.  - 


£<       ar 


&fe+t+4f*e(&>+Z*,]  <34> 


-  *sr*  (<&  (35)- 
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Substituting  Equation  (34)  into  Equation  (35)  gives: 
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Using  Equation  (i.7)  to  eliminate      A3       and  Equation  (23)  to 
eliminate      93        ,   Equation  (36)  becomes  a  relation  involving  only 
one  unknown,       <p3        .     The  solution  must  be  found  by  trial  aid  error. 
In  general,  there  will  be  two  values  of       <p3        which  satisfy  this 
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equation  corresponding  to  the  average  subsonic  and  aupersonic  solutions 


Ttf 
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FIG.  2 -PLOT    OF    MASS     FLOW     VS.    AVERAGE    AXIAL 
VELOCITY    FOR    CHANNEL    FLOW 

If  no  value  of       <P3       satisfies  Equation  (36),  the  prescribed 
mass  flow  cannot  be  forced  through  the  channel  because  of  choking, 
and  the  boundary  configuration  or  upstream  conditions  must  be  altered 
accordingly. 

Having  found  the  average  value  of  the  axial  velocity  at  station  3 , 
an  estimate  of  the  distribution  of  this  velocity  across  the  channel  can 
be  made  considering  the  total  temperature  gradient,  the  blade  shape, 
and  the  pressure  distribution  expected  at  the  downstream  station. 

Using  the  estimated  distribution  of         4>j         ,  Equation  (22)  is 
iterated  across  the  channel  to  determine  the  pressure  gradient.     The 

tangential  velocity  is  found  from  Equation  (17),  the  radial  velocity  from 

eJ.  0 
Equation  (23),  and  the  value  of       .  „  from  Equation  (24).     In  order 

to  evaluate     &3         -^j-^ —       ,  and        ^  ^        at  various  points  across  the 


"/* 
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channel,  it  is  necessary  to  estimate  as  accurately  as  possible  the  slope 
and  curvature  of  each  streamline  at  the  station  in  question.     Knowing 
the  complete  channel  configuration,   slopes  and  curvatures  at  points 
across  the  channel  can  be  easily  estimated  using  a  simple  linear  distri- 
bution scheme,  or  a  large  scale  drawing  can  be  constructed  and  the 
values  of  slope  and  curvature  taken  graphically  from  this.     The  value 
of  ,  ^  in  Equation  (24)  can  be  taken  as  the  average  value  between 

stations  along  the  axis.     However,  it  is  to  be  noted  that  even  a  rela- 
tively large  error  in  the  estimation  of  the  above  quantities  will  have  little 
effect  upon  the  pressure  gradient  obtained  because  of  the  nature  of  the 
integrated  form  of  Equation  (22).     The  actual  pressure  distribution 
across  the  channel  is  found  by  applying  the  energy  equation  at  a  boundary 
to  determine  the  constant  C^. 

With  these  values  of  P*   and        ,  ^  ,  the  axial  velocity 

a.  £3 

integral  can  now  be  iterated  across  the  channel  for  values  of  <>>3  in 
terms  of  the  constant  C ,  •     If  the  curvature  of  the  channel  is  not  large, 

(L03 

the  value  of        .  „  is  very  small  compared  to  the  other  factors  in 

the  integrand  (2G,  for  example),   and  hence  a  relatively  large  error  in 
estimating  has  a  small  effect  on  the  resulting  value  of  4*3* 

The  constant  C  ,    is  obtained  by  iterating  the  mass  flow  relation 
(Equation  (21)  )  across  the  channel,  and  a  new  axial  velocity  distribution 
is  obtained. 

The  above  described  process  is  the  procedure  for  one  complete 
iteration  of  the  given  equations.     Normally  in  the  subsonic  case,  the 
distribution  of  <{>-  thus  obtained  is  used  in  the  next  iteration,   and  the 
process  is  repeated  until  the  axial  velocity  distribution  produced  by  an 
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iteration  is  the  same  as  that  assumed  for  that  particular  iteration. 
This  is  the  solution.     Values  of         T3  and  Q-   are  then  found  by 

use  of  Equation  (13). 

B.     Procedure  for  the  Inverse  Problem 

Generally,  the  same  procedure  of  iteration,  as  described  in 
section  A  above  for  the  direct  problem,   is  used  in  the  solution  of  the 
inverse  problem.     Equations  ( 10),   (13),  (21),  (22),  (24),   (27),  and  (29) 
are  applicable  to  the  isentropic  rotor.     Relations  used  in  the  poly- 
tropic  case  are  derived  in  Appendix  E.     Examples  IX,   X,  a  of 

Appendix  F  are  typical  solutions  to  the  inverse  problem.     One  complete 
iteration  is  oiven  in  Example  IX b  to  illustrate  the  method  used. 

It  i6  noted  that  the  integrand  of  Equation  (29)  contains  a  term 
involving       -jj-         .     Since  in  any  practical  solution  the  axial  velocity 
at  any  point  across  the  channel  will  not  be  zero,  this  term  does  not 
cause  any  real  difficulty  in  the  iteration  procedure.     However,  if  the 
original  estimate  of  the  axial  velocity  profile  at  the  downstream  station 
differs  greatly  from  the  actual  solution,  it  has  been  found  that  succes- 
sive iterations  will  "oscillate"  about  the  real  solution,  but  will 
eventually  converge  as  in  the  direct  problem. 

It  is  obvious  that  the  work  output  function  11  (f3)       ,   given  as 
Equation  (27),  must  be  prescribed  reasonably.     For  any  given  channel 
shape,  blade  speed,   upstream  conditions  and  mass  flow,  there  is  a 
limited  range  of  blade  angles  that  can  be  used  before  the  smooth  flow 
over  the  blades  is  destroyed  by  separation  effects.     The  limitation  thus 
imposed  by  the  maximum  turning  angle  through  the  blade  row  also 
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limits  not  only  the  amount  of  work  that  can  be  extracted,   but  also  the 
distribution  of  this  work  across  the  channel. 

C •     Convergence  of  Subsonic  and  Supersonic  Solutions 

Figure  2  indicates  that,  for  any  given  setup,  there  is  normally 
a  supersonic  as  well  as  a  subsonic  solution  for  the  direct  and  inverse 
problems.     This  is  illustrated  in  Appendix  F  by  Examples  IV  a  and  Va 
for  the  stator  and  Examples  IV  b  and  Vb  for  the  rotor.    In  these 
examples,  for  a  given  boundary  configuration,  blade  spouting  angle, 
and  set  of  upstream  conditions,  both  the  subsonic  and  supersonic 
solutions  have  been  obtained. 

Experience  in  working  various  examples  has  shown  that  if  the 
flow  is  in  the  low  subsonic  range,  point  A  of  Figure  2,  for  example, 
the  iteration  process  converges  very  rapidly,   and  only  two  or  three 
iterations  are  necessary  to  obtain  accuracy  of  four  significant  figures 
in  the  solution.     However  as  the  flow  becomes  highly  subsonic  as  at 
point  B,  the  convergence  becomes  less  rapid,  and  more  iterations  are 
necessary  to  obtain  the  solution.     It  is  interesting  to  note  that,  in  the 
subsonic  range,  the  process  will  cotiverge  from  any  reasonable 
assumed  distribution  of  §*.    It  is  even  possible  to  obtain  a  solution  in 
'  which  the  flow  is  supersonic  in  one  part  of  the  channel  and  subsonic  in 
the  remainder  (Examples  la  and  Ilia). 

It  has  been  found,  however,  that  if  the  flow  is  completely 
supersonic  (point  E  of  Figure  2,  for  example),  the  iteration  process 
diverges  rapidly  away  from  the  solution.     Thus  if  an  axial  velocity 
distribution  is  assumed  whose  average  value  is  represented  by  point  D, 
successive  iterations  will  result  in  decreasing  values  of  ^-i  *nd 
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eventuaily  the  subsonic  solution  at  point  A  will  be  reached.     If  an 
axial  velocity  distribution  corresponding  to  point  F  is  assumed, 
successive  iterations  will  rapidly  increase  93  toward  infinity.     It 
is  the  realization  of  this  phenomenon  which  makes  it  possible  to 
obtain  the  isentrcpic  supersonic  solution  in  a  given  problem.     Having 
made  two  successive  iterations  at  each  of  the  points  D  and  F,  and 
having  noted  the  relative  change  in  the  constant  C,  in  these  iterations, 
a  simple  interpolation  will  give  the  values  of  9,   corresponding  to  the 
solution  at  point  E  once  the  solution  has  been  "straddled". 

It  is  important  to  understand  the  basic  reason  why  the  proposed 
iteration  process  is  successful.    For  either  the  direct  or  inverse 
problems,   six  unknowns  must  be  solved  for  at  the  downstream  station. 
These  are  the  three  components  of  velocity  (    63  ,    X 3  and  9~),  the 
pressure  P*,  the  density  Q-  ,  and  the  temperature      Tj      .     Using 
the  isentropic  (or  polytropic)  relation  along  streamsurfaces,  two  of 
these  (CK   and      T3    )  can  be  eliminated.     The  remaining  four  unknowns 
are  made  functions  of  the  axial  velocity  9,   and  the  boundary  and 
*  streamline  configuration.     The  radial  velocity  is  eliminated  in  terms 
of  9,   by  Equation  (23),  the  tangential  velocity        h3     by  Equation  (17) 
or  (27),  and  the  pressure  P^  by  the  integral,   Equation  (22).     Thus, 
effectively,   only  one  unknown  remains,  the  axial  velocity  9-,  ,  and  for 
any  assumed  distribution  of  9,,   all  other  unknowns  are  determined. 
An  iteration  process  must  be  used  to  determine  9,   because  of  the 
complexity  of  the  equations  involved. 
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D.     Correction  for  the  Displacement  of  Streamlines 

Having  assumed  a  set  of  streamsurfaces  for  any  particular 
example,   and  having  found  the  corresponding  solution  by  the  proposed 
iteration  procedure,  it  is  possible  to  make  a  correction  to  determine 
more  accurately  the  actual  position  of  the  streamsurfaces  and  thus 
refine  the  solution.     The  mass  flow  functions    m     ( /e)       and    m   (£) 
are  plotted  against  radius  for  conditions  at  stations  Z  and  3  as  shown 
in  Figure  6.     Using  this  diagram  and  fixing  the  radial  coordinates  at 
station  3,   it  i3  possible  to  estimate  new  corresponding  upstream 
coordinates.     With  these  corrected  upstream  radial  coordinates,  and 
the  upstream  conditions  corresponding  to  them,  the  problem  is  re- 
worked and  a  new  solution  found.     It  is  to  be  emphasized  that  this  is 
a  small  correction  and  that  the  new  solution  will  not  vary  more  than 
one  or  two  percent  from  the  solution  based  on  the  original  itreara- 
surfaces.     An  illustration  of  this  method  is  given  in  Example  V  b  of 
Appendix  F. 
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IV.     ADDITIONAL  CORRECTIONS  AND  SIMPLIFICATIONS 

A.     Corrections  to  the  Basic  Solution 

If  it  is  desired  to  solve  for  the  flow  pattern  and  distribution  of 
energy  in  the  field  at  any  axial  station  in  the  blade  row  itself,  this  can 
be  done  approximately  by  multiplying  the  axial  velocity  by  a  function 
which  accounts  for  the  blade  thickness  and  thus  reduces  the  effective 
area  of  the  channel.    Assuming  a  blade  thickness  function  t(f)      and 
a  blade  spacing    s  ,  the  increased  axial  velocity  at  any  point  within  the 
blade  row  would  be: 


Mc-H^ro)  <37) 


where    s   and  t  (C)       are  assumed  to  be  known  functions  at  any  radial 
position  across  the  channel.     The  assumptions  used  in  the  direct 
problem,  including  that  of  axial  symmetry,   remain  unchanged,   except 
that  the  blade  spouting  angle  now  becomes  the  mean  camber  line  at 
the  point  in  question.     The  fundamental  equations  and  iteration  pro- 
cedure remain  unchanged  except  for  the  use  of  the  above  indicated 
expression  for  the  axial  velocity. 

One  of  the  basic  assumptions  used  in  the  derivation  of  the 
fundamental  equations  given  in  Part  II,  is  that  the  effects  of  the  vorti- 
city  of  the  flow  downstream  on  upstream  conditions  are  negligible.    A 
correction  to  the  upstream  conditions  due  to  vorticity  effects  down- 
stream could  be  made  using  the  approximate  theory  suggested  by 
Marble  (Reference  i5)  for  the  three-dimensional  case.     If  the  hub 
ratio  is  large,  the  simpler  two-dimensional  approximation  given  by 
Rannie  (Reference  1)  can  be  used.     Once  the  upstream  conditions  have 
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been  changed  using  this  higher  order  correction,  the  corresponding 
new  solution  can  then  be  found  by  the  same  procedure  as  previously 
outlined  in  Part  III. 

®*     Simplification  of  Flow  Patterns 

An  investigation  has  been  made  to  determine  whether  or  not 
there  are  any  special  cases  in  which  the  equations,   used  in  this  analysis 
of  three-dimensional,   compressible  flow,   reduce  to  a  simple  expression 
for  the  axial  velocity.     Consider  the  integrated  radial  momentum 
equation: 


p,  = 


i? 


^(mli-***®** 


■% 


y 

r-i 

(22). 


Since  the  value  of  P-   appears  in  the  integrated  mass  flow 

relation,  Equation  (2i),  to  the  power     -y      ,  it  is  clear  that,   generally, 

no  simplification  can  be  made  for  a  channel  of  arbitrary  shape  where 

the  values  of         G3-T2 —  &n^        &  — rrr-       cannot  be  expressed  in 

df3  r*  <L?3 

terms  of  constants  and  the  radius     ^ 

However,   one  simplification  does  occur  in  the  example  of  a 

rotor  in  a  cylindrical  channel.     In  this  case,  if  the  velocity  after  the 

r 

rotor  is  purely  axial,  Equation  (22)  reduces  to  P,  =  C^-i  .    If,  in 

addition,  the  total  temperature  and  pressure  are  constant  across  the 
channel  at  station  2,  and  if      ^g  =    y~  (vortex  flow).   Equation 

7Z 

(16)  for  a  straight  channel  reduces  to  4?-  =  C,.     The  constant  Cy  *s 
evaluated  by  applying  the  energy  equation  at  a  boundary,  and  the  mass 
flow  integral  becomes: 
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rotor  is  purely  axial,  Equation  (22)  reduces  to  P,  =  C^r-i  .    If,  in 

addition,  the  total  temperature  and  pressure  are  constant  across  the 
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72 
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Ciftf,  (38), 


or,   since     Tzt       is  constant, 


771    = 


*k*    L       Z*  2^j 


"//-tf  09). 


Even  in  this  simplified  case,   Equation  (39)  must  be  solved  by- 
trial  and  error  for  the  constant  C , .     This  particular  example  gives  a 
constant  mass  flow  per  unit  area  and  the  work  is  taken  out  uniformly 
across  the  channel. 

With  the  exception  of  the  above  case,  it  can  be  generally  stated 
that  there  are  no  special  conditions  or  flow  patterns  in  three-dimen- 
sional, compressible  flow  through  turbomachinery  using  this  method 
which  result  in  a  simple  closed -form  solution  for  the  downstream  axial 
velocity. 

C .     Proposed  Optimum  Turbine  Design  Condition 

In  the  design  of  a  turbine,  the  condition  which  normally  controls 
amount  of  heat  that  can  be  released  in  the  combustion  chamber  is  the 
limiting  temperature  that  the  first  stage  rotor  blade  root  can  safely 
withstand  in  continuous  operation.    Since  the  blade  stresses  decrease 
from  root  to  tip,  the  allowable  temperature  is  higher  at  the  tip  than  at 
the  root.     This  fact  allows  the  possibility  of  increasing  the  amount  of 
heat  released  in  the  combustion  chamber,  and  hence  the  energy  in  the 
jet  exhaust,  for  a  fixed  amount  of  work  taken  out  of  a  turbine. 

In  formulating  the  optimum  turbine  design  condition,   reference 
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will  be  made  to  a  one -stage  turbine  with  an  arbitrary,   but  fixed, 
boundary  configuration,   such  as  is  shown  in  Figure  5.     It  is  assumed 
that  the  total  pressure  in  the  combustion  chamber,  the  blade  speed, 
mass  flow  and  total  amount  of  work  to  be  extracted  by  the  turbine  are 
fixed  by  conditions  in  the  compressor.    It  is  further  assumed  that  the 
distribution  of  the  maximum  allowable  total  temperature  across  the 
channel  at  station  2,  based  on  limiting  blade  temperatures,  is  known. 
Since  the  total  enthalpy  is  conserved  through  the  stator  row,  the  limiting 
total  temperature  profile  at  station  1  is  therefore  known. 

In  order  that  the  maximum  energy  be  left  in  the  exhaust  gases 
at  station  3,  the  condition  is  imposed  that  the  velocity  shall  be  axial 
after  the  turbine  (     A^  =  0)  and  that  the  total  enthalpy  shall  be  constant 
across  the  channel  at  this  station  (  **      =  0). 

To  simplify  the  formulation  of  the  proposed  optimum  turbine 
design  condition,  it  is  assumed  that  the  known  limiting  total  temperature 
profile  at  station  1  is  linear,  and  that  the  combustion  chamber  can  be 
constructed  so  as  to  give  this  desired  total  temperature  variation  across 
the  channel. 

The  design  problem  can  then  be  formulated  as  follows:    Given  a 
one -stage  turbine  with  fixed  boundary  configuration,  combustion  chamber 
total  pressure,  blade  speed,  mass  flow,  and  linear  total  temperature 
profile  at  station  I,  what  are  the  blade  spouting  angles  of  stator  and 
rotor  so  that  a  given  amount  of  work  can  be  extracted  by  the  rotor  and 
yet  leave  the  total  temperature  constant  and  the  velocity  axial  at  station 
3? 

With  the  above  assumptions  and  conditions,   Equation  (27)  can  be 
written: 
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"(ft  ~4*K  <40>> 

and  Equation  (10)  in  terms  of  total  temperatures  becomes: 

&»-&.&>*  <4!)- 

If      /_?*-        is  to  be  constant  at  station  3,  then  the  work  must  be  taken 
out  across  the  channel  in  such  a  way  that: 

A  £V  A  %-Tzt  (42). 

Since  the  linear  change  in  total  temperature  at  station  2  is 
assumed  known,  the  distribution  of  work  across  the  channel  is  deter- 
mined by  Equation  (42)  and  the  total  amount  of  work  to  be  extracted. 
The  solution  to  this  simplified  inverse  problem  is  found  in  a  manner 
similar  to  that  illustrated  by  Example  X  of  Appendix  F,   and  the  blade 
trailing  edge  angles  of  stator  and  rotor  are  completely  determined. 

In  view  of  the  discussion  given  in  Part  III,  Section  B,   on 
limiting  blade  angles,  it  will  obviously  not  always  be  possible  to  obtain 
the  above  proposed  optimum  condition,  particularly  if  the  upstream 
total  temperature  gradient  is  large.     Furthermore,  the  above  consi- 
derations could  not  reasonably  apply  to  a  multi-stage  turbine  since 
boundary  layer  effects  and  mixing  would  soon  destroy  the  original 
total  temperature  profile.     However,  the  designer  should  not  overlook 
the  possibility  of  increasing  the  gross  thrust  of  a  turbine  by  the  use  of 
a  total  enthalpy  gradient  based  on  blade  limiting  temperatures.     (See 
Example  XI,  Appendix  F.) 
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V.     RESULTS  AND  CONCLUSIONS 

An  approximate  method  has  been  developed  for  the  solution  of 
three-dimensional,   compressible  flow  in  axial  turbomachines.     The 
methods  and  equations  used  are  extremely  flexible  and  can  be  used  to 
solve  the  direct  and  inverse  problems,   subsonic  or  supersonic,  isen- 
tropic  or  polytropic  for  either  a  turbine  or  a  compressor.     The  upstream 
conditions,   blade  shape  (as  defined  by  the  blade  trailing  edge  angle), 
blade  speed,  and  axial  boundary  configuration  can  be  prescribed  in  the 
direct  problem,  while  the  upstream  conditions,   blade  speed,   blade 
loading,   and  axial  boundary  configuration  can  be  prescribed  in  the 
inverse  problem.     In  each  case,   six  unknowns  are  solved  for  at  the 
downstream  station.     These  are  the  three  components  of  velocity,  the 
pressure,  the  density,  and  the  temperature.     In  the  analysis,   maximum 
utilization  has  been  made  of  the  given  channel  configuration. 

The  equations  are  derived  for  a  blade  row  with  infinitely  many 
blades,   on  the  basis  of  isentropic  (or  polytropic),   non-viscous,   com- 
pressible,  axially- symmetric,   steady  flow.     Radial  or  nearly  radial 
blades  are  assumed  so  that  the  radial  force  can  be  neglected.     Separa- 
tion of  flow  and  boundary  layer  effects  are  not  considered  in  the  analysis. 
For  an  arbitrary  boundary  shape,  a  special  type  of  blading  is  assumed 
whose  only  effect  on  the  flow  is  to  turn  it  through  a  specified  blade 
angle  distribution  across  the  channel.     Using  the  assumption  that  the 
blade  shape  is  specified  by  the  blade  trailing  edge  angle  only,  and 
assuming  that  this  spouting  angle  is  very  nearly  independent  of  condi- 
tions ahead  of  the  blade  row  in  the  operating  range,  the  analysis  of  the 
flow  has  been  made  at  the  station  immediately  downstream  of  the  blade 
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row  in  terms  of  specified  conditions  immediately  upstream  of  the  blade 
row. 

This  method  greatly  simplifies  the  analysis  of  flow  in  turbo - 
machines,   since,  with  the  above  assumptions,  it  is  not  necessary  to 
solve  for  the  complete  three-dimensional  flow  field  through  the  blade 
row,  but  only  to  consider  stations  immediately  upstream  and  down- 
stream of  the  blading. 

The  iteration  process  developed  is  simple  and  straightforward, 
and  the  same  general  procedure  is  used  in  the  solution  of  all  suggested 
types  of  problems.     The  calculations  can  be  done  with  a  slide  rule,   and 
the  accuracy  of  the  solution  can  be  improved  by  increasing  the  number 
of  radial  step3  across  the  channel.     Corrections  can  be  made  for  the 
displacement  of  the  assumed  streamlines  and  for  the  effect  of  down- 
stream vorticity  on  upstream  conditions  to  further  improve  the  solution. 
It  is  believed  that  the  results  thus  obtained  by  this  approximate  method 
are  sufficiently  accurate  for  design  purposes. 

Once  the  boundary  configuration  and  other  necessary  conditions 
have  been  assumed  for  any  particular  problem,   experience  in  working 
various  examples  by  the  proposed  iteration  process  has  shown  that  a 
ba3ic  solution  for  the  direct  problem  accurate  to  four  significant 
figures  can  be  obtained  for  a  stator  in  approximately  four  hours  and 
obtained  for  a  rotor  in  approximately  six  hours.     The  calculations  for 
+>"»  inverse  problem  are  somewhat  simpler  and  less  time  is  normally 

required  for  this  solution. 

oLO 
In  order  to  evaluate  the  term      — —         which  appears  in  both 

<*$ 
the  integral  for  the  axial  velocity  and  the  integral  for  the  pressure  in 
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the  case  of  an  arbitrary  channel  shape,  it  is  necessary  to  estimate  the 
slope  and  curvature  of  the  streamlines  at  the  station  in  question.     It  is 
to  be  noted  that,  for  a  reasonable  channel  shape,  a  relatively  large  error 
in  this  estimation  will  have  little  effect  on  the  solution  obtained.     The 

de 

term  ,  *.  appears  in  the  general  axial  velocity  integral  for  the 

rotor  (Equation  (18)  )  and  for  the  stator  (Equation  (26)  )  as  a  small 
correction  term,   and  a  large  error  in  the  assumed  values  of  slope  and 
curvature  will  have  a  small  effect  on  the  resulting  axial  velocity  profile. 
This  is  especially  true  for  the  rotor  where  the  term  2G  is  normally 
dominant. 

Considering  the  pressure  integral  (Equation  (22)  ),  it  is  to  be 

noted  that  here,   again,  the  contribution  of  the  term      -j-p       is  small 

A2 
for  a  reasonable  boundary  shape.     For  the  stator,  the  term     —j-     is 

the  predominant  one,  and  the  terms        G  -j-j-        and     4>   -jy     are 

usually  small  corrections.     Because  of  the  small  magnitude  of  the  factor 

multiplying  the  integrand  of  Equation  (22)  in  the  iteration,   a  relatively 

large  error  in       —r—        will  have  but  a  slight  effect  on  the  resulting 

pressure  gradient.     After  the  rotor,  where  the  velocity  is  predominantly 

axial,  the  term      —j-~      might  be  of  the  same  order  of  magnitude  as 

Q    ttt7  and      d>  -=— -       .    As  a  result,  it  is  possible  to  have  a 

cLf  T  *? 

negative  pressure  gradient  radially  after  a  rotor  if  the  slope  and  curva- 
ture of  the  channel  is  large. 

It  is  impossible  to  make  any  general  statement  as  to  the  relative 
contribution  of  the  integral  (Equation  (18)  )  and  the  constant  of  integration 
C.   to  the  resulting  axial  velocity  profile.     The  relative  size  of  the 
integral  and  the  constant  may  be  of  any  magnitude  depending  on  the 
channel  shape,  total  temperature  gradient,   blade  spouting  angles,   and 
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upstream  conditions  of  any  particular  example. 

Several  interesting  observations  are  possible  from  the  results 
of  the  examples  given  in  appendix  F.     The  channel  shapes  assumed  for 
these  iterated  proolems  are  given  as  Figures  3,   4,  and  5.     In  each 
example,  the  assumptions,   channel  configuration,   upstream  conditions 
radial  coordinates,  and  other  data  are  specified,   and  the~resulting 
solution  at  the  downstream  station  is  tabulated.     One  complete  iteration 
is  given  in  Examples  Vb,  VIII  b,   and  IX  b  to  show  the  method  used  in 
each  case. 

Examples  I  and  II  use  the  channel  configuration  shown  in 
Figure  3.     In  both  of  these  examples  of  the  direct  problem,  the  mass 
flow,  blade  speed,   blade  spouting  angles,  and  total  heat  content  in  the 
flow  are  the  same.     Example  I,  however,  has  no  radial  total  tempera- 
ture gradient  at  station  1 ,  while  Example  II  has  a  large  radial  enthalpy 
gradient.     The  effect  of  a  total  temperature  gradient  across  the  channel 
in  changing  the  axial  velocity  profile  is  clearly  evident  from  the  results 
obtained.     Normally  for  a  stator,  it  is  to  be  expected  that  the  positive 
pressure  gradient  would  make  the  axial  velocity  higher  at  the  root  than 
at  the  tip.     This  is  the  case  in  Example  la  for  no  temperature  gradient. 
However,   Example  II a  shows  that  the  effect  of  the  positive  total 
gradient  imposed  is  to  make  the  axial  velocity  higher  at  the  tip  than  at 
the  root,   and  reduce  the  xvlach  number  slightly  at  all  stations  across  the 
channel.     Also,  the  resulting  total  pressure  at  station  3  has  been  in- 
creased by  the  use  of  the  radial  enthalpy  gradient. 

Examples  IV  and  V  use  the  channel  configuration  of  Figure  5, 
and  have  the  same  mass  flow,   blade  speed,  blade  spouting  angles,  and 
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conditions  at  station  1.     Example  IVa  is  the  subsonic  stator  solution 
while  Example  Va  is  the  supersonic  solution  for  the  same  stator. 
Examples  IVb  and  Vb  are  the  subsonic  and  supersonic  rotor  solutions, 
respectively,  for  conditions  after  the  subsonic  stator  of  Example  IVa. 
The  mass  flow  functions  for  the  subsonic  and  supersonic  rotor  examples 
are  plotted  in  Figure  6.    From  this  graph,  the  position  of  the  displaced 
streamlines  has  been  estimated  and  plotted  on  Figure  7.     This  plot  of 
the  streamline  shift  in  each  case,   clearly  indicates  that  the  shift  is  in 
one  direction  when  the  "governing"  Mach  number  is  subsonic,   and  in  the 
opposite  direction  when  the  "governing"  Mach  number  is  supersonic. 
The     governing"  Mach  number  has  been  defined  by  Monroe  (Reference  7, 
page  51)  as  the  Mach  number  based  on  the  meridional  velocity  instead 
of  the  total  velocity.     The  results  shown  in  Figure  7  clearly  agree  with 
Monroe's  observation  that  the  deflection  of  the  " meridional  Mach 
surfaces"  is  in  different  directions  for  subsonic  and  supersonic 
"governing"  Mach  numbers,  and  that  the  deflection  of  the  streamsurfaces 
will  be  less  if  the  governing  velocity  is  transonic  than  if  it  is  entirely 
subsonic  or  entirely  supersonic. 

Examples  VI,  VII,  and  VIII  use  the  channel  configuration  of 
Figure  5,  the  same  conditions  at  station  1,   and  the  same  blade  speed, 
blade  spouting  angles,   and  mass  flow.     Example  VI  is  the  isentropic 
case,  Example  VII  is  a  polytropic  solution  with    n  =  1.  37  (77  =  94.  5  per- 
cent), and  Example  VIII  is  a  polytropic  solution  with    n  =  1.  33 
(  ~rj     -  86.  8  percent).     From  these  examples,  the  effect  of  a  decrease 
in  "polytropic  efficiency"  on  the  three-dimensional  velocity  field  and 
distribution  of  energy  in  the  flow  can  clearly  be  seen.     When  the 
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"efficiency"  has  dropped  to  94.  5  percent,  the  average  axial  velocity  at 
station  3  has  increased  2.  7  percent,  the  average  total  pressure  has 
decreased  2.  5  percent,  and  the  work  output  has  increased  4.  7  percent 
when  compared  with  the  isentropic  solution.     As  the  "polytropic 
efficiency"  drops  to  86.  8  percent,  the  results  of  Example  VIII  show  that 
the  average  axial  velocity  at  station  3  has  increased  8.  2  percent,  the 
average  total  pressure  has  decreased  8.  0  percent,   and  the  work  output 
has  increased  17.0  percent  when  compared  with  the  isentropic  solution. 

The  results  of  these  examples  clearly  show  the  reason  for  the 
"loading  up"  of  the  first  stages  of  a  multi-stage  turbine  as  the  efficiency 
drops.     As  the  efficiency  decreases,   more  work  is  taken  out  in  the 
initial  stage,   and  less  total  pressure  is  left  for  the  remaining  turbine 
stages.     Thus,  the  last  stage  of  a  multi- stage  turbine  may  well  act  as 
a  compressor  if  the  efficiency  of  the  system  has  dropped  sufficiently. 
In  view  of  this  phenomenon,  it  is  obvious  that  the  design  engineer  should 
design  his  machine  for  the  expected  efficiency  in  the  operating  range, 
and  not  on  the  basis  of  isentropic  flow. 

Example  XI  has  been  constructed  to  illustrate  the  increase  in 
gross  thrust  obtainable  by  the  use  of  a  radial  total  temperature  gradient. 
Assuming  a  one -stage  turbine  configuration  similar  to  that  shown  in 
Figure  3,  and  keeping  the  mass  flow,  combustion  chamber  total  pressure, 
blade  speed,  and  work  extracted  constant  in  each  case,  the  gross  thrust 
of  the  jet  exhaust  has  been  calculated  with  and  without  a  total  tempera- 
ture gradient.     In  order  to  give  the  example  more  physical  meaning, 
blades  were  used  which  were  tapered  in  cross-sectional  area  from  hub 
to  tip.     Having  calculated  the  centrifugal  stress  at  various  radial  points 
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along  the  blades,  and  using  the  high  temperature  data  for  "Vitallium" 
(U.S.  21)  (NB-10)    given  in  References  16  and  17,  it  -was  estimated 
that  the  rotor  blade  root  could  safely  stand  a  temperature  of  1450  F, 
while  the  blade  tip  temperature  could  be  1850  F. 

Using  this  data,  the  gross  thrust  of  the  turbine  was  calculated 
for:    (a)  A  constant  radial  total  temperature  of  1450  F,  and  (b)  a  linear 
total  temperature  gradient  of  400  F  with  the  same  root  total  tempera- 
ture of  1450  F.     The  results  of  these  examples  show  that  an  increase 
in  gross  thrust  of  917.8  pounds,   or  7.7  percent,  was  obtained  by  the 
use  of  the  radial  total  temperature  gradient  indicated. 
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APPENDIX  A 

SYMBOLS  AND  DEFINITIONS 

8&  -  Stator  spouting  angle  measured  from  the  plane  perpendicular 

to  the  axis. 

^3  -  Rotor  spouting  angle  measured  relative  to  the  rotating  wheel. 

c  -  Axial  distance  between  stations  in  feet. 

C,  -  Constant  obtained  from  axial  velocity  integral. 

Cy  -  Constant  obtained  from  radial  momentum  integral. 

C  -  Specific  heat  at  constant  pressure. 

P 

A  -  Used  to  denote  incremental  changes  of  a  quantity  in  the  radial 

direction. 

p(|2)      -  Cotangent  of  stator  spouting  angle  (3a  as  a  function  of  radius. 

G  (fs)      "  Cotangent  of  rotor  relative   spouting  angle    f33    as  a  function 
of  radius .        G(L)  =    *3  +  f* 

Y  -  Ratio  of  specific  heats. 

g  -  Acceleration  of  gravity  (32. 17  ft/sec   ). 

H(yj)      -  Dimensionless  worn  output  of  rotor  as  a  function  of  radius. 

Woo 
m  -  Dimensionless  mass  flow.       m  =  2Trgp  /?Q 

§ 
M  -  Mach  number  based  on  total  velocity.    M  =  v. 1 

Mg  -  "Governing"  Mach  number  as  defined  on  p.   87  of  Reference  7. 

M3-\ 


4>^+e&' 


n  -  Polytropic  exponent. 

n  -  Polytropic  efficiency.  [         JV   -  n  _ 

to  -  AnguJar  velocity  of  rotor  in  radians  per  second. 

p  -  Pressure  in  pounds  per  square  foot. 

P  -  Dimensionless  pressure.   P  -   -^r 


p 

Q 
R 


2         4 
Density  (lb.  sec.     /ft.    ). 


czRZf> 


r 

* 

T 

T 

u 

e 

V 

w 

♦ 

I 

W 
z 

f 


3*-f 


-  Dimensionless  density.        Q  ~       -p 

ft2 

-  Gas  constant  (1715 ). 

2o_, 
sec     F 

•  Reference  radius  in  feet. 

-  Radius  from  the  axis  in  feet. 

-  Dimensionless  radius,  f  -   ~^~~ 

-  Absolute  temperature  in  degrees  Rankine. 

-  Dimensionless  temperature.     ?   -    ^a  /^e  T 

-  Radial  component  of  absolute  velocity  in  feet  per  second. 

-  Dimensionless  radial  component  of  absolute  velocity.    0=  ^  ^ 

-  Tangential  component  of  absolute  velocity  in  feet  per  second. 

-  Dimensionless  tangential  component  of  absolute  velocity.    A  =  ^  % 

-  Axial  component  of  absolute  velocity  in  feet  per  second. 

*v 


Dimensionless  axial  component  of  absolute  velocity.    <j>  -  ^  g 

Dimensionless  total  absolute  velocity.   §  =   \  ^  +  A  +  & 

Mass  flow  in  pounds  per  second. 

Axial  distance  along  axis  of  rotor  from  reference  point  in  feet. 

Dimensionless  axial  distance.  £  -  -jh- 

Subscript  to  denote  station  immediately  upstream  of  stator. 

Subscript  to  denote  station  immediately  downstream  of  stator 
and  immediately  upstream  of  rotor. 

Subscript  to  denote  station  immediately  after  rotor. 

Subscript  to  denote  inner  boundary  of  channel. 

Subscript  to  denote  outer  boundary  of  channel. 

Slope  of  streamline. 

Rate  of  change  of  slope  of  streamline. 
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APPENDIX  B 
Note:    In  Appendices  B,  C,  D,   £  and  F,  Equations  of  lower  number 
than  43  refer  to  relations  previously  given  in  the  text.     The 
numbering  of  equations  is  consistent  throughout  the  thesis. 
Development  of  Equations  for  Direct  Problem  (Isentropic). 

The  equations  are  derived  for  a  rotor  with  infinitely  many  blades 
on  the  basis  of  isentropic,  non-viscous,  compressible,  axially-symmet- 
ric,   steady  flow.   Radial  or  nearly  radial  blades  are  assumed  so  that 
the  radial  force  can  be  neglected.     The  channel  shape,  biade  shape, 
mass  flow,  and  upstream  conditions  are  assumed  specified  in  any  par- 
ticular example. 

The  derivation  follows  the  general  procedure  given  in  Refer- 
ence 1,  Section  2:2.  2  for  the  incompressible  case.   The  energy  equa- 
tion is  written  for  any  two  points  (for  example,  b2  and  b~,  Fig.    5)  on 
a  streamline  surface  between  stations  2  and  3.     Using  absolute  veloc- 
ities, the  equation  in  dimensionless  form  is: 


■&rn-4.**-£r*,*T  k *  ( V AA *»' 


-  T  $s  (43) 


where 
r 


-r-i 

and 


tn-    £rt*±(6l+£+  ?t)t 


The  energy  equation  is  now  differentiated  in  respect  to  the  rad- 
ius at  stations  2  and  3: 
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Using  the  gas  equation,     L3  -  -x-    ,  the  first  term  of  the  right 
side  of  Equation  (44)  becomes,  with  the  iscntropic  relation: 


■y-,  df,         r-i  <tf,\  Q3j        *-'  d(3  [  <pa    r3     I 


Using  the  differential  form  of  the  continuity  equation  for  a 
stream  cylinder  of  small  radial  extent, 

and  differentiating  Equation  (45)  as  a  product,  the  following  is  obtained: 

.2L  ±-  (&.)   -   J_  d8  .  _2L  _*    A_  (  f& )  <Ps  (s  n« 

r-t  dfs  \  Q3j  Q3  dfs  +  r-i    p^  <t£  \<?a  J  <Pef&  UD'- 

Using  the  radial  momentum  equation  at  station  3, 

and  combining  Equations  (8),  (11),  (15)  and  (44),  the  following  differ- 
ential equation  results: 
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Qs 


hdf?     *'  P2rdf2\  tel<Pe4e  *  *  cLf3W+*3  +  VsJ  +  cL  fs 

\r-,    *$,     dfe     J*& 


=  (Jk?4l&  \jr_  dT^_     *(fe\)     _r   fe  %_    a_  fp/\ 


dfs 


(46). 


For  the  direct  problem  it  is  assumed  that  the  rotor  spouting 
angle    fij     is  known  as  a  function  of  radius  at  station  3.     It  is  there- 
fore possible  to  describe  the  relation  between    A3  and     <f>a     along 
the  radius  with  a  function  G{f3)  such  that: 


^-*%M-& 


(17). 


Eliminating  Aj  from  Equation  (46)  with  this  relation  the  fol 
lowing  differential  equation  is  obtained: 

—       r  ■"         -0 

_  (PsY '  4>3$3        r  dTat      dtfahe)       r  r  (p3\r  d     In    P/)\    .  X  9  &- 

\Pa.J  3I&  1  *t  <?£,     «*£      ~  r-i  c4pz)   dfe  I*  <?a  J\    f*d  £ 
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G% 


_fPsf  fa  j"   r  dTtt    atfM      r  r(Pa)*9  (o  f£\ . .  9&.  p/9 


9<Ps 


+  <P3 


Psf     fa 


PJ    **h 


y  an*  _  d(f*K) 


-  £,QpJ  afj^Qj 


>+^i-aG    (47). 
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If  P,,  fjy    -f^3     an<*  the  known  quantities  at  station  2  are 
considered  to  be  functions  of  $3  (to  be  determined  in  the  iteration 
process)  the  integral  for   <f>3    can  be  immediately  written  down  as: 


4, 


A" 


*        *^ 


I       (  PX   6 


-r-t 


£7z\ii)  dj&i;) 


*«*(-*#; 


fs*t,  +  C, 


(13). 


An  alternate  form  for  the  square  bracket  in  the  above  inte- 
grand can  be  obtained  by  using  the  energy  and  radial  momentum  equa- 
tions at  station  2: 


'%-A%+i#*<i-^-^(ff-i)i&$[  <i9» 


In  the  case  that  P,.  «  1  (isentropic  expansion  to  station  2  from 
stagnation  conditions),  the  bracketed  term  in  the  integrand  may  be 
further  simplified  to: 


£3ft(>-*V-£i-(W 


The  integral  form  of  the  mass  flow  equation  is: 

7i 


(20). 


(21). 


and  the  integral  form  of  the  radial  momentum  equation  is: 
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r-i 
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Equations  (10),   (13),  (17),  (18),  (21),  (22),  and  (24)  are  used 
in  the  solution  of  the  direct  problem  for  the  rotor  (isentropic). 
In  the  case  of  a  stator,   Equation  (17)  becomes: 


**  =  nfj*. 


(25), 


G  is  replaced  by  F  in  Equation  (13)  and  the  terms 


and 


2G  are  omitted.     Equation  (18)  becomes  Equation  (26)  for  the  isentropic 


stator.     Equations  (13),   (21),  (22),  and  (24)  (with  appropriate  sub- 
scripts) remain  unchanged,  and  the  energy  relation  becomes  Equation 
(9). 
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APPENDIX  C 
Development  of  Equations  for  Direct  Problem  (Polytropic). 

The  equations  for  a  polytropic  process  through  a  rotor  are 
derived  with  the  same  assumptions  as  used  in  Appendix  B  with  the 
exception  that  the  isentropic  relation  (13)  is  now  replaced  by: 

t  ■  m-  at 

where  n  is  related  to  the  so-called  polytropic  efficiency  77    for  an 
expansion  (n  <  -y)  by:    (see  Reference  12,  p.  449) 

-2L.77   =     -21  (49). 

72-/   /  Y-l 

The  energy  equation  (10),  the  continuity  equation  (8),  the  rad- 
ial momentum  equation  (11),  and  the  blade  relation  (17)  remain  the 
same  in  differential  form. 

The  first  term  of  the  right  side  of  Equation  (44)  now  becomes: 

*•'  HAQsl        *-i    Qsdfs  *-<   <?/  dfs 

Y  d_  (_fi* \   =    J_  dft    ,   _/ /_  d% r    Ps  dQs 

r-,  *£(<?,/         Qj  dfs        r-i    Q3  cL$3       r-t  Qj  <tfs  . 


But, 

_i i_dR?       r    P*    dQs    =       r  p*?  d.     I  Pjr\ 

t--/    Q3d.f3        r-i  q*  df3  r->3     *$3  \   Q3  I. 

So,  as  before, 

-y  a    (P?\.    j_  cL3      _y_p¥<L-  (#l)  /so). 
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The  right  side  of  the  differential  equation  is  obtained  in  the 
same  form  as  Equation  (46)  with  a  correction  term  as  follows: 


Pi 


Qs 


■¥)» 


n-r 
rn 


_ZL n'r 'sL-  (£l)  _    JLD? 


dfs 


■n-r 
m. 


(51). 


Using  the  differential  form  of  the  continuity  equation  (8),   Equation 
(51)  finally  becomes: 


72 


(52). 


The  differential  equation  corresponding  to  Equation  (46)  for 
the  isentropic  case  now  is: 


_  JLp^JL-  f^jl  +  cfi  ?*>+     ™    (6?) 


d  P3™ 


(53), 


*$3 


The  relation  corresponding  to  Equation  (47)  becomes: 


+  4 


72-/ 


The  integral  equation  corresponding  to  Equation  (18)  is: 
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-Q*  \Pel 


4s 


<fiaf* 


(55). 


If  P,t  =  1  at  station  2,  the  square  bracket  under  the  integral 


will  be: 


(56). 


The  integrated  form  of  the  mass  flow  equation  is 


772     = 


% 


and  the  integral  form  of  the  radial  momentum  equation  is: 

4. 


R  = 


"^  l^!/  fa   ^  *>*<3)disc' 

—  7i 


_  _?2_ 
72-/ 


(57) 


(58). 


Equations  (10),  (17),  (24),  (48),   (55),  (57)  and  (58)  are  used 
in  the  solution  of  the  polytropic  direct  problem  for  the  rotor. 

In  the  case  of  the  stator,  Equation  (25)  is  used  in  place  of 
Equation  (17),  and  the  terms      Z~f(^z  ^z)    am*  2G  are  omitted  in 
Equation  (55).     Equations  (10),  (24),  (48),  (57)  and  (58)  remain  un- 
changed (with  appropriate  subscripts). 
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APPENDIX  D 
Development  of  Equations  for  the  Inverse  Problem  (Isentropic). 

The  following  equations  are  derived  with  the  same  assump- 
tions as  given  in  Appendix  B  with  the  exception  that  now  the  blade 
relation.  Equation  (17),  is  replaced  by  the  following  given  function 
which  describes  the  work  output  of  the  rotor  as  a  function  of  the 
radius   j3    : 


H(f3)  =  t*K+t,K 


(27). 


For  the  isentropic  case,  the  derivation  remains  unchanged 
up  to  an  including  Equation  (46),  which  may  be  written  as: 


i*  ihfa*>kH*)  -  &M 


4,     *  z>4, 


98 


r.p*¥2-(3L 


*-i  3       ^  I  <?x 


(28). 


Using  Equation  (27)  to  eliminate    A3,  the  differential  equa- 
tion corresponding  to  (47)  is  developed  as  follows: 
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(59) 
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(60). 


Assuming  that  the  right  side  of  the  above  Equation  (60)  is 
a  known  function  of   £3      (to  be  determined  in  the  iteration  process), 
the  integral  equation  for  the  axial  velocity  can  be  written  as: 


*- 


(PsY  4s  J    r   dT2t        r  rfyd  (  p/\     (H-t*h*)<L(t  A  ^ 


d.G3 


'  <P3  ('+    4/   U43 


df3+c, 


(29). 


Since  the  mass  flow  integral  and  the  radial  momentum  inte- 
gral remain  the  same.  Equations  (10),  (13),  (21),  (22),  (24),  (27)  and 
(29)  are  used  in  the  solution  of  the  isentropic  inverse  problem  for 
the  rotor. 
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APPENDIX  E 
Development  of  Equations  for  the  Inverse  Problem  (Polytropic). 

The  folio-wing  derivation  is  made  assuming  the  same  conditions 
for  a  rotor  as  given  in  Appendix  C .     The  work  output  of  the  rotor  as  a 
function  of  radius  is  again  given  by  Equation  (2?). 

The  differential  equation  (53)  still  holds,  which  is: 

r-i  *     *&  \<?2  l\v'd?3  +(r-0(n-o  \Qjd4p     J  (53) 


The  resulting  integral  equation  is  similar  to  Equation  (29)  with 


the  addition  of  a  correction  Jterm: 
.6 
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(r-,H~-,AQzlcL$3  ('*     f/    /df 
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The  mas 8  flow  and  radial  momentum  equations  remain  the 
same  as  in  Appendix  C,  and  Equations  (10),  (24),  (27),  (48),  (57), 
(58)  and  (61)  are  used  to  solve  the  polytropic  inverse  problem  for 
the  rotor. 
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APPENDIX  F 
Iterated  Examples 

The  following  data  are  used  for  Examples  I-X  inclusive: 


u 

Y 
R 


=  734  rad/sec  (7000  RPK.)      pu  =  16,300 


lb 


ft' 


R     =  1.2500  ft. 
o 


=  1.40 


=  1715 


ft" 

"To", 
sec      * 


u  R     =917.5  ft/sec 
o 


o>2R2   =   .8418  x  106  ~ 


C     =   .2396 


BTU 

lb°F 


c  =  .250  ft 


sec 


Example  la . 

Solution  for  the  stator  of  Fig.  3  (cylindrical  channel),  with  no 
radial  total  temperature  distribution  and  coordinates  for  equal  areas 
in  the  plane  perpendicular  to  the  axis  .     This  is  a  solution  for  the  direct 
problem  (F  ( ^a  )  given)  for  the  isentropic ,  subsonic  case. 

For  the  stator  with  F  constant,  P..  =  I,  and  9  =  0,  Equations  (18), 
(21),  (22),  (13),  (25)  and  (43)  of  Appendix  B  become: 


A- 


i'^O+F*) 


6    ± 


J4i 


(ffil^^O-^i^H 


6 


171     = 


i 


*,[%)' *kf**(* 


+  ^hs    (62) 


(63) 


P^ 


r~i      i 


y    ?> 


J    f* 


<±fe    +  fl 


r 

7-1 


(( 4) 


i-iff-tir 


(65) 
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K  =  Ffr 


(F  =  const.) 


(66) 


r-i  L'* 


£r£+i(&*f) 


(67) 


M   = 


§ 


W 


(K/Iach  number) 


(68). 


Conditions  at  station  I:    (given) 


Pu=1.00 


dTn 


It 


■// 


F 


=  0 


.9723 
F  =  2.75 


It 


m  = 


23828 

Z/t=  4.1968             <J>j  =  .4857 

23354 

T,  =  4.1631              \,  =    0 

.020417 

W  =  114.7  -i£-    T14.  =  206 

Solution  at  station  2: 


"2 

6 

.8000 

+2 

.7926 

* 

2.1796 

r2 

.4928 

Zi 

3.4283 

C2 

.14372 

M 

1.059 

Ar2 

+.0001 

.6508 

C,  = 

.8169 

u2 
.8544 

.7479 
2.0567 

.5365 
3.5128 

.15274 

.986 
-.0002* 


TO  -     P 

*2t      Mt 


"2 
.9055 

.7105 
i.9539 

.5730 
3.5794 

.16007 

.929 
-.0001 


sec 


2t 


d2 
.9539 

.6785 
1.8659 

.6041 
3.6340 

.16623 

.880 
-.0004 

a  C 
W  It 


0 


c2 
1.0000 

.6508 
1.7897 

.6307 
3.6739 

.17144 

.841 
-.0002 

2=0 
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Note:      The  column  headed  aT.  above  represents  the  difference 
between      T  _  as  calculated  by  the  iteration  process  and  T^ 
obtained  by  applying  the  energy  equation  along  the  stream- 
lines indicated  after  the  solution  has  been  reached.    The  same 
applies  in  the  solution  of  the  following  examples  .     This  is  a 
rough  check  on  the  accuracy  of  the  solution.    Obviously,  if  an 
exact  solution  had  been  found,  A  T        would  be  zero. 

Example  lb. 

Solution  for  the  rotor  of  Fig.  3  (cylindrical  channel)  with  no 

radial  total  temperature  gradient  and  radial  coordinates  for  equal  areas  . 

This  is  an  isentropic,  subsonic  solution  for  the  direct  problem  (G    {73    ) 

given)  for  a  rotor  . 

For  the  rotor  with  G  constant  and  P-,.  =  1  ,  the  applicable 

relations  are  Equations  (63),  (64),  (65)  (with  appropriate  subscripts) 

and  the  following  three  equations: 


*,m 


6* 


0+G*) 


f*[*w-#-4*v 


dn 


+86 


*  4i  (69) 

if* 


^J  =  e%~^    <G  "  const.) 


(70) 


■y-t  "*r      '2  "&   ~    r-l"3'     Z    l   "J       rM  J       *73'l3  I71)- 

Conditions  at  station  2  are  the  same  as  after  the  stator ,  Example  la. 


G  =  1.411 


d.Tzt 


* 


=  0 


m  =  .020417 
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Solution  at  station 

3: 

a3 

b3 

c3 

d3 

c3 

£ 

.8000 

.8544 

.9055 

.9539 

1.0000 

4>3 

.6624 

.6779 

.6955 

.7146 

.  7346 

*3 

.  1346 

.  1021 

.0759 

.0544 

.0365 

P3 

.5862 

.5865 

.5865 

.5865 

.5865 

?3 

3.6025 

3.6034 

3.6034 

3.6035 

3.6035 

Q3 

.16269 

.16278 

. 16276 

. 16276 

.16277 

1  3t 

.6244 

.6254 

.6264 

.6273 

.6284 

*W 

3.6679 

3.6699 

3.6717 

3.6734 

3.6751 

Q3t 

. 17019 

. 17041 

. 17061 

.17077 

.  17099 

bX3 

0 

Cx  «  .  5710 

-.0007 

-.0017 
C2  =  .8585 

-.0035 

-.0058 

Example  II a. 

Solution  for  the  stator  of  Fig.  3  (cylindrical  channel)  with  the 
same  conditions  and  assumptions  as  given  in  Problem  la  except  with  a 
radial  total  temperature  gradient. 

Equations  (62)  -  (68)  of  Example  la  are  applicable. 


Conditions  at  station  1:    (given) 


I 


x, 


it 


Q 


It 
It 


.8000 

3.3319 
1.00 

.29569 

.9723 

3. 3548 
.28981 

.4355 

0 

F  =  2.75 


Ml 

.8544 

3.8252 

1.00 
.26142 
.9723 

3. 7946 
. 25622 

.4628 

0 


"I 

.  9055 

4.2415 

1.00 
.23577 
.9723 

4.2076 
.23108 

.4871 

0 


m  a  .020417 


ul 

.9539 

4.6361 

1.00 
.21570 
.9723 

4. 5990 
.21141 
.5096 
0 
cLT/t 


cl 
1.000 

5.0118 

1.00 
.19953 
.9723 

4.9717 

.19556 

.5298 

0 


H 


=  +8.  1495 
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Soluti 

on  at  station 

2: 

a2 

fit 

b2 

C2 

d2 

e2 

& 

.8000 

.8544 

.9055 

.9539 

1.0000 

<t>z 

.6778 

.6805 

.6829 

.6848 

.6864 

A* 

1.8639 

1.8714 

1.8780 

1.3832 

1.8876 

P2 

.5295 

.5713 

.6061 

.6355 

.6611 

?a 

2.8200 

3.2599 

3.6762 

4.0729 

4. 4527 

Q2 

. 18777 

.17525 

. 16488 

« 

.15602 

.  14845 

M 

.999 

.930 

.883 

.842 

.805 

A?^ 

0 

-.0012 

-.0052 

-.0104 

-.0172 

C 

.  =  . 5565 

C2  =  .8339 

P2t  *  Plt  ' 

r   =  r 

•Q2t  = 

Qlt  ' 

82=0 

Example  lib 

Solution  for  the  rotor  of  Fig.  3  (cylindrical  channel)  with  the 
•came  conditions  and  assumptions  given  in  Problem  lb  except  with  a 
radial  total  temperature  gradient.     Equations   (63) ,  (64),  (65),  (69) »  (70) 
and  (71)  of  Example  I  (with  appropriate  subscripts)  are  applicable.     Condi" 
tions  at  station  2  are  the  same  as  after  the  stator,  Example  Ha,  above. 


G  =  1.411, 
Solution  at  station  3: 

*& 

-=  3.1495, 

m  =  .020417 

a3 

b3 

c3 

d3 

e3 

is 

.8000 

.8544 

.9055 

.9539 

1.0000 

*3 

.5377 

.6174 

.6920 

.7628 

.8303 

h 

-.0413 

.0168 

.0709 

.1224 

.1716 

P3 

.6009 

.6009 

.6009 

.6009 

.6012 

?3 

2.9237 

3.3072 

3.6670 

4.0081 

4.3334 

Q3 

.20551 

.18168 

.16386 

. 14990 

.13871 

P3t 

.6313 

4 

.6381 

.6422 

.6447 

.6461 

Q3t 

.21289 

.18964 

.17184 

.15759 

.14603 

r3t 

2. 9653 

3. 3643 

3.7373 

4.0895 

4.4235 

tx3 

0 

•f.0035 

+.0012 

-.0039 

-.0126 

Cl 

= 

.4635 

C2  =  .8646 
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After  correcting  the  upstream  conditions  at  station  2  for  the 
displacement  of  the  streamlines  from  the  position  originally  assumed, 
the  following  corrected  solution  at  station  3  is  obtained:    (See  Problem 
Vb  for  method.  ) 


<3 
*3 
*3 

'-3 

r3 

°3 
P3t 

t3t 

°3t 

at, 


tt3 

8000 

.  5311 
-. 0506 
.  6029 
2.9268 
.20602 
6327 
2.9675 
.21325 
0 


u3 
8544 

.6106 

0072 

6029 

3.2771 

. 18411 

.6378 

3.3303 

19166 

-.0001 

Cj  =  .4578 


c3 
.9055 

.6853 
0615 

.  6029 
3.6296 

. 16623 

.6417 
3.6951 

. 17381 
-. 0022 


C2  = 


u3 
9539 

7560 

1128 

.6029 

3. 9868 

. 15127 

.6453 

4  0648 

15879 

-  0055 

8654 


"3 

1.0000 

.8234 
.  1618 
.6032 

4.3374 
. 13917 
.6476 

4.4263 
. 14641 

-.0118 


Example  Ilia- 

Solution  for  the  stator  of  Fig.    4  (divergent  channel  with  straight 
walls)  with  a  radial  total  temperature  gradient.     This  is  an  isentropic, 
subsonic  solution  for  the  direct  problem  for  a  stator  (F  (£2  )        given). 

For  the  stator  with  F  constant  and  P.     =  1,   Equations     (63)  ,  (65) 
and  (66)  are  applicable  in  addition  to  the  following  four  relations: 


A  = 


Fz 


£'~W; 


4 


[p>'  ti  I/-'  *ft  (     z  '1   *&. 


7+F' 


frn 


0 


& 


ELE  {12.) 


t    61 


«- 


£ 


r-i 


(73) 


-  £r.*i(fr£  +  4t) 


(74) 


3& 


it  m  *  40i 


(75), 


Conditions  at  station  1:    (given) 


d, 


f, 

.80 

.85 

.90 

.95 

1.00 

r,t 

3.9931 

4.0950 

4. 1968 

4.2987 

4. 4006 

Qu 

. 25043 

. 24420   • 

.23828 

.23263 

.22724 

pu 

1.00 

1.00 

1.00 

1.00 

1.00 

pi 

.9696 

.9696 

.9696 

.9696 

.9696 

T, 

3.9531 

4.0591 

4.  1600 

4.2610 

4. 3620 

Ql 

. 24497 

.23888 

.23309 

.22756 

.22229 

*i 

.4944 

.5014 

.5069 

.5112 

.5142 

9i 

-.0247 

.0000 

.0253 

.0511 

.0771 

*i 

0 

0 

0 

0 

0 

F-2.75  ^-^     =+2.0375        m  =.  024639       W  =  138.  4  — 

dl  sec 


c  =  .  25  ft 
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Solution  at  station  2: 


u 

0a 
P2 

Q2 

M 

a* 


.790 

.7410 

2.0377 

37C 

.  5247 
3.3212 

. 15798 
1.004 
- 


Cx  =  .6017 


"2 

850 

.7046 
1  9376 
0 
.  5716 
3.4904 
16378 
.930 
-.0027 


Cz  =  .8317 


"2 

•  910 

.6735 
1    8521 

.  0337 

.612?. 
3.6479 

.  16732 

.870 
-    3061 


"2 

•  970 

.6469 
1.7790 

.0647 

.6472 
3.7966 

. 17051 

.817 
-. 0103 


c2 

1.  30 

.6236 
1. 7149 

.0935 

.6777 
3.^376 

. 17210 

.779 
-.0138 


P2t=Plf   °2t  =  Qlt,    ^  =  ^ 


Example  Illb. 

Solution  for  the  rotor  of  Fig.  4  (divergent  channel  with  straight 
walla)  with  a  radial  total  temperature  gradient.  This  is  an  isentropic, 
subsonic  solution  for  the  direct  problem  for  the  rotor  (  G    (  ^3)     given). 

Equations    (63),  (65),  (7C) ,  (73) ,  (75)  (with  appropriate  sub- 
scripts) and  the  following  two  relations  are  applicable: 

4 


l"°U*) 


irH 


+  ■%   (76) 


£«*  -k  V  £  5  *  t  (*>  K*  4H,  4. 


(77). 


Conditions  at  station  2  are  the  same  as  after  the  stator,   Example 
Ilia  above. 


G  ■  1.411       m  =  .024639 


dn* 


*4* 


-  +  1.6979       c  =  .25  ft. 
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Solution  at  station  3: 


h 

P3 

Q3 

?3 

p 


3t 
Q3t 

&T3 


*3 

.780 

.  4567 

-.0228 

-. 1356 

.6500 

. 18409 

3.5308 

.6712 

. 18836 

3-5634 

+.0001 


u3 
.850 

.4987 
0 
-. 1463 

.6500 

.  17954 
3.6209 
6749 

. 18442 
3.6600 
+.0005 


Cl* 


3871 


^3 

.920 

.5423 
.0271 

-1548 
.  650? 
17525 

3.7117 
.6780 
. 18051 

3.7560 

-. 0013 


C2S 


"3 

990 

.  5873 
.0587 

-'.  1613 
.6508 
. 17119 

3. 8027 
.6804 
. 17672 

3.8514 

-  0048 

8842 


c3 

1 .  060 

.6347 
.0952 

-.  1644 
.6508 
. 16720 

3.8923 
.6828 
. 17304 

3.9459 

-.0092 


Example  iVa 

Solution  for  the  stator  of  Fig.    5  (divergent  curved  channel  walls) 
with  a  radial  total  temperature  gradient.     This  is  an  isentropic,   sub- 
sonic solution  for  the  direct  problem  for  a  stator  (F  (^z)  given).     The 
same  equations  as  used  in  Example  Ilia  are  applicable. 


Conditions  at  station  1:    (given) 


t, 


It 


It 


Q 


It 


*1 
ti 

Pi 

♦l 
°1 
h 


F  =  2 


*1 

.80 

3.7894 

1.00 
. 26389 
.9400 

3.7231 
. 25248 
.6804 

-.0340 
0 

75 


ul 

.84 

3.9524 
1.00 
.25301 
.9400 
3.8832 
. 24207 
.6960 
0 
0 
dTit 


di 


.88 

4. 1152 

1.00 
. 24299 
.9400 

4.0433 
. 23248 
.7090 
.0354 
0 

=  4. 0746 
c  =  . 


"1 
.92 

4.2783 

1.00 
.23374 
.9400 

4. 2034 
.22363 
.7205 
.0720 
0 

m  =  . 029242 

250  ft. 


cl 
.96 

4.4413 

I.  00 
.22516 
.9400 

4.3636 
.21542 
.7293 
.  1094 
0 


Al 

1.00 

4. 6043 

1.00 
.21719 
.9400 

4. 5237 
. 20780 
.7365 
.  1473 
0 


y;   =  164.2  — 
sec 
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Solution  at  station  2: 

a2 

b2 

c2 

d2 

e2 

f2 

f* 

.785 

.842 

.899 

.956 

1.013 

1.070 

+* 

.7051 

.6763 

.6531 

.6344 

.6192 

.6072 

*z 

1.9390 

1.8598            1 

.7960 

1.7446 

1.7028 

1.6698 

92 

-.0705 

.0135 

.0914 

.1649 

.2353 

.3036 

P2 

.5417 

.5860 

.6222 

.6518 

.6761 

.6962 

r« 

3.1805 

3.3928            3 

.5937 

3.7860 

3.9713 

4.1519 

o 

V  2 

.17032 

.17272 

.17313 

.17217 

.17022 

.16767 

M 

.977 

.909 

.854 

.809 

.776 

.747 

blz 

0 

0 

.0013 

-.0039 

-.0069 

-.0118 

Cj  =  .5694 

C2  =  .8393 

*2t  =    Plt    *    ^        Tl±      '  C2t  =       It 


Example  IVb 

Solution  for  the  rotor  of  Fig.  5  (divergent  curved  channel  walls) 
with  a  radial  total  temperature  gradient.     This  is  an  isentropic , 
subsonic   solution  of  the  direct  problem  for  the  rotor  (G  (&    )      given)  . 
The  same  equations  as  used  in  Example  Hlb  are  applicable. 

Conditions  at  station  2  are  the  same  as  after  the  stator  , 
•Example  IVa . 


cLT 
G  =  1.730  -T^-     =2.8593  m  =   .029242 
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Solution  at  station  3: 


6 

+3 

G3 
P3 


lJr 


Q3t 

P3t 

M 


*3 

.760 

.  1958 
-.4213 
-.0294 

.6948 
3.4149 

.  20346 
3.4660 

.21117 

.7320 

.212 
0 


u3 

.848 
.2641 

-.3911 
.0106 
.6981 

3.5669 
.  19573 

3. 5998 
.  20027 
.7209 
.210 

+.0010 


Cl  = 


c3 

.936 
.3260 

-.3720 
.0750 
.7010 

3.7180 
. 18852 

3.7535 
. 19303 
.7246 
.219 

-.0003 

1594 


d3 
1.024 

.3852 
-.3576 

.  1617 

.7027 
3.8682 

. 18167 
3.9064 

. 18619 

.7274 

.236 
-.0049 

=  .9012 


c3 
1.  112 

.4458 

-.3408 
.2719 
.7010 

4.0126 
. 17466 

4.0568 
. 17950 
.7283 
266 

-.0113 


1.200 
.5117 

-.3148 
.  4094 
.6927 

4. 1461 
16707 

4.2018 
. 17272 
.7257 
.302 

-.0196 


Example  Va. 

Supersonic  isentropic  solution  for  the  stator  of  Example  IVa 
(divergent  curved  channel  walls  and  radial  total  temperature  gradient) 
The  conditions  at  station  1,   the  data,  and  the  applicable  equations  are 
the  same  as  for  Example  IVa. 

Solution  at  station  2: 

c 


+2 

#a 
P2 

Q2 
M 

tX9 


~2 

•  9303 

2.5583 
-.0930 

.3171 
2.7294 

. 11617 
1.367 
0 


Cx  « .7512 


"2 
.8906 

2.4491 
.0178 
.3715 

2.9784 
. 12471 

1.274 

+.0039 


C2  = J203 


2 

.8601 

2.3653 
.  1204 
.4169 

3  2051 
. 13007 

1.  188 

+.  0034 


P      *  P 
*2t 


"2 

.8371 

2  3020 
.2176 
.4547 

3.4157 
. 13310 

1.  123 

-.0012 


v2 
.8201 

2.2553 

.3116 

4856 

3.6131 
.13440 

1.080 

-. 0032 


*2 

<  8087 

2.  2239 

.4043 

.5102 

3.7990 

. 13430 
1.038 
-.0178 


It' 


2t      ^ 


It' 


t   =  r 


/t 


66 


Example  Vb. 

Supersonic  isentropic  solution  for  the  rotor  of  Example  IVb 
(divergent  curved  channel  walls  and  radial  total  temperature  gradient)  . 
The  conditions  at  station  2  are  assumed  the  same  as  after  the  subsonic 
stator  of  Example  IVa .     The  data  and  equations  used  are  the  same  as  that 
used  for  Example  IVb. 

One  complete  approximation  for  this  problem  is  given  below  to 
show  the  method  of  iteration.    The  procedure  of  correcting  for  the  dis- 
placement of  the  assumed  streamlines  is  then  given  to  show  the  method 
for  this  correction. 

The  iteration  is  started  by  assuming  an  axial  velocity  distribu- 
tion across  the  channel.     The  tangential  velocity  is  then  calculated  using 
Equation  (70).     The  value  of         ,    J  is  found  from  Equation  (75)  and 

the  pressure  distribution  from  Equation  (73)  and  the  application  of  Equa- 
tion (77)  at  the  inner  channel  boundary.     The  axial  velocity  distribution 
is  then  determined  using  Equation  (76)  and  the  constant  is  determined  by 
the  mass  flow  Equation  (63) .    The  iteration  is  continued  until  the  values 
of  <J>,  assumed  agree  with  those  resulting  from  the  iteration  process.     It 
is  important  to  note  that  in  the  supersonic  case,  the  iteration  diverges 
rapidly  away  from  the  solution.    Thus  if  an  axial  velocity  distribution  is 
assumed  that  is  greater  than  the  solution,  successive  iterations  will  in- 
crease $,  until  the  pressure  F-  has  reached  zero.     If  an  axial  velocity 
distribution  is  assumed  which  is  less  than  the  solution,  successive 
iterations  will  result  in  decreasing  values  of  §-  and  eventually  the  subsonic 
solution  will  be  reached.     The  supersonic  solution  is  found  by  a  simple 
interpolation  once  the  solution  has  been  "straddled'  . 
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The  reverse  ia  true  for  the  subsonic  solution  for  a  given  setup, 
The  iteration  process  will  converge  toward  the  solution  from  any 
reasonable  assumed  distribution  of  $,. 

Data  used  in  the  solution: 

*    cLTa 

,2 


&f3    =.088  -Zj  Jj^*«"  10. 0075        G=  1.730         2G  =  3.460 


=  .7496 


r? 


a3       b3      c3      d3      e3      f3 
^  .760      .848    .936    1.024     1.112   1.200 


(if.) 

Factor  .02117     .02125  .02133   .02138   .02144 


-.150      .040    .230     .420     .610    .800 
-.250      .100     .450     .800     1.150   1.500 


1.4310     1.5415  1.6411   1.7305   1.8098 


L—  1.2283     1.1316   1.0509    .9824    .9235   .8722 


t-7496 

*2     .2 


Q  (fi-fi-,)  .01299     .01357  .01489   .01694   .01719 

^2    2 


-fyfa) 


-.7702     -.8526   -.9333  -1.0018   -1.0842 


y-^r- (fi^^-t)        .006498    .006235  .005993  .005768  .005560 


%t 


0± 
("Factor"  above  =  — ^ /  £  '+Gedf 

Try  C.-  1.3920      ^"^% 
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a3  b3  c3  d3  e3  f3 

<p3  1.7098  1.6727  1.6943  1.7698  1.9004  2.0912 

G<f>3  2.9580  2.8938  2.9311  3.0618  3.2877  3. 617S 

X3  2.1980  2.0458  1.9951  2.0378  2.1757  2.4178 

e3  -.2565  .0669  .3897  .7433  1.1592  1.6730 

X3  4.8312  4.1853  3.9804  4.1526  4.7337  5.8458 

-^-  6.3568  4.9354  4.2525  4.0553  4.2569  4.8715 


/, 


5.6461         4.5939        4.1539        4.1561         4.5642 


6 

~63  -.0948    .2283    .5665    .9512    1.4161 


3.6750    3.6681    4.0181    4.7261    5.8336 


-Q?h.  f-,3484        -.8374      -2.2763      -4.4955      -8.2680 


~J'J3 


d<t>. 


*-         5.0235        4.9820         5.2060         5.6770        6.4060        7.4200 


9f, 

..7535    .1993    1.1974    2.3843    3.9077    5.9360 


bf3\*fI  ' 

m3  - 


j  -  -j 

4274  .1673  .7624  1.4158  2.1855  3.1368 

§%r       -1.1809  .3666  1.9598  3.8001  6.0932  9.0728 

-<k¥r  +2.0191  -.6132  -3.3205  -6.7254  -11.5795  -18.9730 


~4sW  f.7030         -1.9668       -5.0229      -9.1524      -15.2762 


r-i 


&%  *  +.04352       -i-   .01116     -    .01885   -   .05475    -      .10553 


69 


The  constant  C^  is  found  by  applying  Equation  (77)  at  the  inner 

channel  boundary: 
T3  =*  3.  7894-.  14286  [4.  8312+2.  9234+.  0658+2(1.  5221+1.  6705)]  =  1.  7600 
P3  =  .541 7(.  5534)7^7   *  .  0683  C2  =  (.  0683)^=  .4645 


r-/ 
pr 

Ps 

m 

Average 

[l-P3  7  ) 
Average 


a3 
.4645 

.0683 

.  1261 

.2278 

.5355 


D3 
5080 

0934 

1594 

2694 


c3 
5192 

1008 

1620 

2725 


u3 
5003 

0886 

1359 

2404 


e3 
4456 

0591 

0874 

1754 


x3 
3400 

0229 

0329 

0373 


2486         .  2709 

.4920  .4309 

5137         .4864 


.2564  .2079  .1313 

.4996  .5543  .6601 

.4902  .5269  .6072 


The  axial  velocity  integral  is  iterated  in  steps  across  the  channel  as 
follows:  4^  ^ 

.02117  [.3557(5.  1414-. 7702)+3.  46-.  4073]  =  .09754  .09754 
.02125  [.4176(4.  8681  -.8526)  \*.  46  +  1.1633]  =  .13388  .23142 
.02133  [.  4208(4.9062-.  9333)  f3.  46  (-2.  8302]  =  .17090  .40232 
.  02138  [.  3598(5.  2735-1.  0013)+3.  46+4.  947o]  =*  .21260  .61492 
.02144  [.2376(6.  0766-1.  0842)+3. 46+7.  583^=  .26220  .87712 
The  mass  flow  equation  is  iterated  to  solve  for  C,  as  follows: 
.  029242=.  003229(.  04877+1. 1799CX)+.  003676(.  16448+1.  0912Cj) 
+.  003813(.  31687+1.  0166Cj)+.  0Q3522(.  50862+.  95290^ 

^.  002257(.  74602+.  8978Cj) 

C     -  .023790  .  ,    ,Q?0 
Cl      7517*58?"  1,392° 


Average 
04877 

.16448 
.31687 
. 50862 
. 74602 


Solution  at  station  3: 
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4. 
e3 

P3 

T3 

Q3 

Zit 

°3t 

P3t 

M 

M 

g 

aT3 


a3 

1.7098 
2. 1980 
-.2565 

.0683 
1.7601 

. 03880 
2.8773 

. 13256 

.3815 
1.730 
1.100 
K0001 


°3 
1.6727 
2.  0458 

.0669 

.0934 
2. 0079 

. 04653 
3.0093 

. 12794 

.3849 
1.570 

.993 
K  0028 


^3 

1.6943 

1.9951 
.3897 
.  1008 

2.1365 
.04718 

3. 1204 
. 12162 
.3795 

1.540 

1.005 

-.0170 


a3 

1.7698 

2.0378 
.7433 
.0886 

2. 1406 
.04139 

3.2056 
. 11358 
.3641 

1.620 

1.110 

-.0548 


c3 

1.9004 
2.  1757 
1.1592 

.0591 
1.9793 

. 02986 
3.2572 

. 10373 

.3379 
1.880 
1.340 
-. 1065 


*3 

2.0912 
2.4178 
1.6730 

.  0229 
1.5653 

.01464 
3. 2649 

.09198 

.3000 
2.420 
1.810 
-.1603 


Correction  of  supersonic  solution  for  displacement  of  streamlines: 

4 


772 


7£ 


Kits) 


°2 

785 
0 


u2 

.842 
. 005498 


"2 

.899 
.011201 


m 


ts 


<ti  ~-  f3Q3<P3$3^ 


*3 

760 


b3 

.848 


w3 

.936 


w2 

.956 
.017071 


°3 

1.024 


c2       *2 
1.013     1.070 
.023095   .029242 


e3 
1.112 


A3 

1.200 


.005258       .011609       .018374         .024639       .029242 


The  functions  m(^)  and  m($3  )  are  plotted  against  radius  on 

Fig.   6.     From  this  graph,  a  corrected  ^  is  obtained  and,  hence, 
new  conditions  at  station  2. 
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Corrected  values  at  station  2: 


& 

<Pa 
<9a 

r« 

P2 


*2 

.785 

.7051 
1.9390 
-.0705 
3.7394 

.5417 
3.1805 

.17032 


°2 

.339 

.6778 
1 . 8640 

.0090 
3.9438 

.5837 
3.3815 

. 17259 


u2 

.902 

.6322 
1.7926 

.0955 
4.1240 

.6241 
3. 6043 

.17315 


Corrected  solution  at  station  3: 


*3 
h3 

83 

P3 

r3 

Tst 

°3t 

P3t 

a7^ 


*3 
1.7398 
2.  2499 
-.2610 

.0605 
1.7003 

. 03558 
2.8660 

. 13104 

.3643 
+.0006 


Cx  =  1.4165 


u3 

1.7007 

2.0942 
.  1020 
.0842 

1.9450 
.04330 

2.9896 
.12665 
.3790 

+.0033 


C2  =  .4487 


u3 

1.7151 

2.0311 
.2573 
.0953 

2. 1067 
. 04523 

3.1139 
.12076 
.3763 

-.0069 


u2 

.968 

.6305 
1.7338 

.1804 
4.3127 

.6580 
3.8266 

.17197 


°3 

1.7686 

2.0357 

.4422 

.0959 

2.2072 

.04344 
3.2376 
.11325 
.3666 
-.0368 


c2 

1.026 

.6157 
1.6933 

.2514 
4.4785 

.6816 
4.0135 

.16978 


c3 

1.8722 

2. 1269 
.9361 
.0733 

2.1223 
.03452 

3. 3064 
. 10449 
.3458 

-.0881 


*2 

1.070 

.6072 
1.6698 

.3036 
4. 6043 

.6962 
4. 1519 

.  16767 


x3 

2.0493 
2. 3453 
1.6394 

.0295 
1.6832 

.01754 
3. 2898 

.09375 

.3071 
-.1032 


Note:  Examples  VI,   VII  and  VIII  use  the  channel  configuration  of  Fig. 
5  and  the  same  conditions  at  station  1.     Example  VI  is  the  isen- 
tropic  solution  (y  ■  1.  40),   Example  VII  is  a  polytropic  solution 
for  n  a  1.  37,  and  Example  VIII  is  for  n  =  1.  33.     From  these 
examples  it  is  possible  to  see  clearly  the  effect  of  a  decrease 
in  polytropic  efficiency  on  the  solution  for  the  same  mass  flow, 
upstream  conditions  and  channel  configuration. 
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Sxample  Via. 

Solution  for  the  stator  of  Fig.    5  (divergent  curved  channel 
walls)  with  a  radial  total  temperature  gradient.    This  is  an  isentropic 
(v  «  1.  40)   subsonic  solution  for  the  direct  problem  for  a  stator 
(F(^2)  given).     The  radial  coordinates  have  been  chosen  for  equal 
areas  in  the  plane  perpendicular  to  the  axis.   The  same  equations 
as  used  in  Example  Ilia  are  applicable. 


Conditions  at  s 

tation  1: 

(given) 

al 

bl 

cl 

dl 

el 

fl 

i 

.8000 

.8438 

.8854 

.9252 

.9633 

1.0000 

?,t 

3.7894 

3.9679 

4.1374 

4.2995 

4. 4543 

4. 6043 

Plt 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Qlt 

. 26389 

.25202 

.24170 

. 23253 

. 22443 

.21719 

Pl 

.9500 

.9500 

.9500 

.9500 

.9500 

.9500 

°1 

. 25439 

.24295 

.23300 

. 22421 

.21640 

. 20937 

r, 

3.7343 

3.9102 

4.0772 

4.2369 

4. 3900 

4.5373 

4>t 

.6202 

.6355 

.6481 

.6583 

.6657 

.6716 

9, 

-.0310 

.0030 

.0367 

.0701 

.1026 

.1343 

F  ■  2.75          m 

=  .  026865       %%*- 

=  4. 0746 

\  = 

0 

Solution  at  station  2:    (v  ■  1. 40) 


i 


z 
<Pz 

9a 

P2 
°2 

M2 
&T2 


a2 
.7850 
.5772 
1.5873 
-.0577 
.6712 
. 19848 
3.3814 
.775 
0 

Lat  L/1r 


w2 
.8496 
.5513 

1.5161 
.0193 
.7095 
.19723 

3.5974 
.754 

-.0013 

P   =  P 
^2t 


"2 

.9098 
.5317 

1.4622 
.0865 
.7378 
.19451 

3.7930 
.675 

-.0025 

If  °2t  s 


1. 


3. 


It 


"2 
9661 
5166 
4206 
1453 
7592 
19103 
9738 
66S 
0037 
C 


c2 

1.0194 
.5046 

1.3876 
.1986 
.7763 
.18732 

4. 1437 
.630 

-.0061 

=  . 4661 


*2 
1.0700 

.4948 
1.3607 

.  2474 

.7898 

. 13349 
4.3041 

.609 
-.0080 
C2  =  .8923 
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Example  VIb 

Solution  for  the  rotor  of  Fig.   5  (divergent  curved  channel  walls) 
with  a  radial  total  temperature  gradient.     This  is  an  isentropic  (7=1.40) 
subsonic  solution  for  the  direct  problem  for  the  rotor  {G(f  )  given). 
The  radial  coordinates  are  for  equal  areas.   The  same  equations  as 
used  in  Example  Illb  are  applicable. 

Conditions  at  station  2  are  the  same  as  after  the  stator,  Example 
Via. 


G  =  1.96  m  =  .026865 

Solution  at  station  3:    (v  =  1.40) 


d?i 


2t 


'£ 


=  2.8593 


4 


3 
<p3 

o3 

P3 

r3 

°3 

*n 

P3t 

°3t 

ttf3 


*3 
.7600 
.1652 

-. 0248 

-.4362 
.7550 

3.4970 
.21589 

3.5281 
.7788 
. 22073 

-.0001 


w3 
.8661 
.2370 
.0187 

-.4016 
.7595 

3.6683 
. 20707 

3.6993 
.7824 
.21150 

-. 0002 


v3 

.9605 
.2982 
.0844 

-.3760 
.7626 

3.8290 
.19916 

3.8605 
.7847 
.20326 

-.0024 


a3 

1 .  0464 
.3543 
.1660 

-.3520 
.7635 

3.9802 
.19181 

4.0126 
.7853 
.19572 

-.0072 


c3 

1.1258 
.4089 
.2616 

-. 3244 
.7606 

4.1197 
. 18460 

4.1549 
.7837 
. 18859 

-.0135 


*3 

1.2000 
.4658 
.3726 

-.2870 
.7532 

4. 2460 
.17738 

4.  2867 
.7788 
.18167 

-.0219 


Cj  =  .1329 


C2  ■  .9228 


Example  Vila 

Solution  for  the  same  channel  configuration,  mass  flow,  and 
conditions  at  station  1  as  for  Example  Via,  except  for  a  polytropic 
process  (n  =  1.  37  ,     V  =  94.  5  /o)  through  the  stator. 

Equations  (43),  (57),  (58),  (66), (74),  (75)  and  the  following 
integral  are  applicable  for  the  polytropic  stator  (F  ■  const): 


<k- 


74 


±~JWm^('-n*i 


I    '!■•'•}, 


a  02 

r2 


(78). 


(r-7L)     £ 


72-/ 


/£ 


(r-o (■».-/)    <p£  dfz 

F  =  2.75  m«  .026865 

Solution  at  station  2:    (n  =  1.  37) 

b 


c 


P2 


2t 


Q 


2t 

*r2 


"2 

.7850 

.5897 

1.6216 

-.0590 

.6450 

3.3635 

.19176 

.9789 

. 25834 

+.0001 


'2 

.8496 

.5652 

1.5543 

.0203 

.6837 

3.5778 

.19110 

.9821 

.24751 

-.0017 


"2 

.9098 

.5469 

1 . 5040 

.0890 

.7127 

3.7726 

.13889 

.9845 

.23793 

-.0022 


u2 
.9661 

.5330 

1.4657 

.  1499 

.7347 

3.9526 

.18585 

.9863 

.22938 

-.0038 


-El 


=  4.0746 


c2 
1.0194 

.5222 

1.4360 

.2055 

.7523 

4.1218 

.18251 

.9874 

.22164 

-.0066 


*2 
1.0700 

.5135 

1.4121 

.2567 

.7662 

4.2814 

.17895 

.9882 

.21461 

-.0090 


Cx  =  . 4762 


C2  »  .8883 


Example  Vllb 

Solution  for  the  same  channel  configuration,  blade  shape, 
and  mass  flow  as  for  Example  VIb,  except  for  a  polytropic  process 
through  the  rotor  (n  ■  1.  37,  J    «  94.  5°/o). 

Conditions  at  station  2  are  the  same  as  after  the  stator, 
Example  Vila. 

Equations  (48),  (57),  (58),  (70),  (75),  (77)  and  the  follow- 
ing integral  are  applicable  for  the  polytropic  rotor  (G  =  const, 
Pzt*  l>: 


t- 


75 


r  A 


i'H-o') 


(Pel  6  /    *■*   <*f         *-*  %    *£  ^  fo  f    <tV'*  e) 


*fjf'   *% 


*&  '    <P3  (f-Ofa-'AQeldf 


i    dfs 


,   c, 


+8G 
(79). 


G  a  1.96        ma  .026865 
Solution  at  station  3:    (n  =  1.  37) 


a  2.8593 


h 

P3 

r, 

Q3 

*W 

P3t 

°St 

&T9 


*3 
.7600 
.1664 

-.4339 

-.025a 
.7387 

3.4890 
.21172 

3.5199 
.7620 
.21646 
0 


°3 
.8661 
.2496 

-.  3769 
.0197 
.7430 

3.6590 
. 20306 

3.6839 
.7608 
.20653 

-. 0044 


"3 
.9605 
.3120 

-. 3490 
.0883 
.7457 

3.8190 
. 19524 

3.8422 
.7617 
.  19823 

-.0092 


u3 
1.0464 

.3637 
-.3335 

,1704 

.7463 
3.9694 

.18799 
3. 9947 

.7632 

.19102 
-.0137 


c3 
1.1258 

.4097 
-.3228 

.  2621 

.7439 
4. 1094 

.18101 
4. 1404 

.7636 

.18443 
-.0177 


*3 

1.2000 
.4514 

-.3153 
.3611 
.7379 

4. 2382 
.17410 

4.2807 
.7640 
.17849 

-.0194 


Cx  a  .1338 


C2  =  .9215 


It  is  interesting  to  note  that,  compared  to  the  isentropic  case 
(Examples  VIa,b),  the  polytropic  flow  through  the  identical  stator 
and  rotor  results  in  the  following  percentage  changes  at  station  3 
(maximum); 


$3:  5. 3°/o 
^3:  f7.7°/o 
P3:    -2.3°/o 


T3:  -0.3% 
P3t :  -3.  0% 
Work  out:  +6.6°/o 
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Sxamplc  Villa. 

Solution  for  the  same  channel  configuration,  mass  flow,  and 
conditions  at  station  1  as  for  Example  Via,  except  for  a  polytropic 
process  (n  =  1.33,  V  =  86.3°/o)  through  the  stator.     The  equations 
used  in  Example  Vila  are  applicable. 


F 

■  2.75 

m  ■  .026865 

3746 

Solution  at  station  2: 

(n  =  1.33) 

a2 

b2 

c2 

d2 

e2 

*2 

U 

.7850 

.8496 

.9098 

.9661 

1.0194 

1.0700 

4>t 

.6225 

.5994 

.5828 

.5704 

.5613 

.5545 

** 

1.7119 

1.6483 

1.6027 

1.5686 

1.5436 

1 .  5249 

6a 

-.0622 

.0216 

.0948 

.1605 

.2209 

.2772 

P2 

.5378 

.6280 

.6584 

.6817 

.7000 

.7146 

r2 

3.3149 

3. 5286 

3.7225 

3.9022 

4.0695 

4. 2279 

°2 

.17731 

.17796 

.17685 

. 17470 

.17199 

.  16902 

P2t 

.9387 

.9470 

.9532 

.9571 

.9608 

.9631 

P2t 

. 24770 

.23863 

.23035 

.22262 

.21564 

.20918 

&T£ 

-.0001 

-.or  02 

-.0019 

-.0044 

-.0071 

-.0107 

5,= 

"  r*     cx  - 

.5027 

C2  =  .8765 

Example 

vnib 

Solution  for  the  same  channel  configuration,  blade  shape,  and 
mass  flow  as  for  Example  VIb,  except  for  a  polytropic  process 
through  the  rotor  (n  «  1.33,  7  ■  86.8  Jo). 

Conditions  at  station  2  are  the  same  as  after  the  stator,  Ex- 
ample Villa.     The  equations  used  in  Example  Vllb  are  applicable. 
One  complete  approximation  for  the  polytropic,  subsonic,  direct  prob- 
lem for  the  rotor  is  given  below  to  demonstrate  the  iteration  procedure, 


Data  used  in  the  solution: 
m  =  .026865  G  =  1.96 


rr-T^"  10.0075 

r-i  *f 
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2G=  3.9200  .,  =  .7935 

1+G^ 


m 

brv3 

Factor 


A  6 


.7600     .8661     .9505    1.0464    1.1258    1.2000 


-.1500     .0791     .2829     .4684     .6398     .8000 


-.2500     .1720     .5474     .8891    1.2049    1.5000 


.02083     .01874    .01716    .01593    .01495 


1.6273     1.7561    1.8551    1.9330    1.9949 


9.0490    8.9897    8.9174    8.8593 


9.6600     9.8311   10.0646   10.2042   10.2714 


(r.-n.)    (P£\      2.0529     2.1450    2.2328    2.3169    2.3975 
(Y-/)(n-j)\Qz  J 

1.2433    1.1208    1.0325     .9646     .9102     .8653 


793  f 


& 

pffrf;-')  .01532  .01530         .01516         .01495         .01470 

2^/_^  W.£  )      .006801  .006099       .005064       .004509       .004072 

(For  definition  of  'Factor"  above,  see  Problem  Vb) . 
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Try  Cj  =  .1568 

a3                  b3                  C3  d3  C3  f3 

4>3                .1949           .2717  .3312  .3775  .4164  .4524 

G4>3              .3820            .5325  .6492  .7399  .8161  .8867 

/Vj             -.3780  -.3336  -.3113  -.3065  -.3097  -.3133 

93             -.0292     .0215  .0937  .1768  .2664  .3619 


A 

7, 


1880     .1285     .1009     .0898     .0852     .0818 
.1582      .1147     .0953     .0875     .0835 


is 

&3  -.0039  .0576  .1352  .2216  .3141 

-GL^T-  +.0019  -.0418        -.1308        -.2501         -.4043 

*73 

ih/jf.)    +.3207        -.1296        -.3559        -.4518        -.4635        -.4084 
9*3  [*t)3 

d>(£-£.\     -.0487     .0467     .1813     .3356     .5017     .6786 


983 
'** 

903 


+.2720         -.0829         -.1746         -.1162         +.0382         +.2702 


+.0946  -.1287        -.1454        -.0390        +.1542 


*  "'  -.022  +.0388        +.0  +  .015  -.0670 


9*3 

***** 

Arj   *■  ^.00094         +.00068      +.00008      -.00066       -.00158 
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The  constant  C,  is  found  by  applying  Equation  (77)  at  the  inner 
channel  boundary. 

^3  =  3.  7894-.  14286  [.  1429 +.  0380  ►.  0009+2(1.  3438-.  2873)1   =  3.4616 


n-1 


P3»  .  5878(1.  0443)1 


7000      C2  *(.7000)    n    =  .9153 


l3 


n-1 


P3  n         .9153  .9162  .9169  .9170  .9163  .9148 


.7000  .7028 

1.1909         1.1191 


.7049  .7052  .7031 

1.0706        1.0345        1.0044 


*3 

<P^ 
^2 

P     I 
(^)n     1.1404         1.0883         1.0526        1.0258         1.0033 

Av. 


.6984 
.9773 

.9829 


n-1 
n 


'     a.      n  * 


1.1143         1.0704         1.0392         1.0145        .9931 
.9157  .9165  .9169  .9166        .9155 

+.0380         +.0227         +.0026         -.0209      -.0490 


The  axial  velocity  integral  is  iterated  in  steps  across  the  channel  as 

follows:  ,  A 

^  4$        <P3       Average 

.02083[l.  8139(9. 1313-3.  8461)+3.92f.  0946^.07861  a  .09602  .09602  -04801 

13267 

.  01874[l.  8797(9.  0490-9.  0107)+3.  92-.  1287  (-.04871=:  .07331  .16933  ' 

19781 

.01716[l.  9278(3.  9397-9.  2287)^3. 92-.  1454^00581  =  .05697  .22630* 

r  4  .25001 

,01593[l. 9610(3. 9174-9. 3537)+3. 92-. 0390-.0484|=  .04742  .27372 

29523 

.  01495[l.  981 1(8.  8593-9. 4040) +3.  92+.  1 542-.  1175]  *  .04302  .31674  ' 


The  mass  flow  equation  is  iterated  to  solve  for  C,  as  follows: 


30 


.  026865=.  017071(.  04801  -H.  1820CX)  f.  016377(.  13267  +1.  0766C}) 
f.  01 5754(.  19781 +.99850^+.  015167(.  25001 +.93740^) 
+.  014597(.  29523+.  8877^) 


wl  "  .GS6717       "w 

3 

Solution  at  station  3:    (n  =  1. 

33) 

a* 

b3 

c3 

d3 

c3 

f3 

& 

.7600 

.8661 

.9605 

1.0464 

1.1258 

1.2000 

4>3 

.1949 

.2T17 

.3312 

.3775 

.4164 

.4524 

h 

-.3780 

-.3336 

-.3113 

.3065 

-.3097 

-.3133 

e, 

-.0292 

.0215 

.0937 

.1768 

.2664 

.3619 

p3 

.7000 

.7028 

.7049 

.7052 

.7031 

.6984 

t} 

3.4617 

3.6285 

3.7862 

3.9354 

4.  0740 

4.2038 

C3 

. 20220 

. 19367 

.18615 

.17921 

.  17275 

.16518 

T3t 

3.4876 

3.6503 

3.8062 

3.9582 

4.1049 

4. 2456 

P3t 

.7185 

.7177 

.7181 

.7204 

.7219 

.7229 

Q3t 

. 20600 

. 19659 

.18863 

.18183 

.17605 

. 16929 

t>Z3 

-.0001 

-.0047 

-.0108 

-.0155 

-.0178 

-.0202 

Compared  to  the  isentropic  case  (Examples  VIa,b),  the  poly* 
tropic  flow  with  n  ■  1.  33  through  the  same  stator  and  rotor  results  in 
the  following  maximum  percentage  changes  at  station  3: 


<J>3:  +18,0% 
X^\  +20.8% 
P3:  -  8.3% 


r3:  -1.1% 

P3t:  -9.3% 
Work  out:  +  19.5% 


Example  IXa 

Solution  for  the  stator  of  Fig.  5  (divergent  curved  channel 
walls)  with  a  radial  total  temperature  gradient  and  radial  coordinates 
for  equal  areas.   This  is  an  isentropic,  subsonic  solution  for  the  direct 
problem  for  the  stator.     The  resulting  solution  at  station  2  is  used  in 


8i 


Example  IXb  as  the  upstream  conditions  for  a  solution  of  the  inverse 
problem  for  a  rotor. 

The  equations  used  in  Example  Ilia  are  applicable. 

Conditions  at  station  1   :     given 


/, 


It 


It 


Q 


It 


r, 

e, 
*, 


i 

.8000 
3.7894 
1.00 

. 26389 

.9400 
3.7231 

. 25249 
-.0340 

.6804 


wl 
.8438 

3.9679 

1.00 
.25202 
.9400 

3.8985 
.24113 
.0033 
.6970 


cl 
.8854 

4. 1374 

1.00 
.24170 
.9400 

4.0650 
.23126 
.0403 
.7103 


F  ^  2.75 

Solution  at  station  2: 


m  a  .029242 


ul 
.9252 

4.2995 

1.00 
.23258 
.9400 

4. 2243 
. 22253 
.0768 
.7215 

dTit 


cl 
.9633 

4.4548 

1.00 
. 22448 
.9400 

4. 3768 
.21478 
.1125 
.7303 


f, 


=  4.  0746 


1.0000 

4. 6043 

1.00 
.21719 
.  9400 

4. 5237 
.20781 
.1473 
.7365 


P2 

°z 


a2 

.7850 

.7008 
1.9272 
-.0701 

.5460 
3.1379 

.17126 
0 


T  -  Z 

L3t  ~  L/t 


Cx  a  . 5659 


w2 
.8496 
.6688 

1.8392 
.0241 
.5954 

3.4217 
. 17402 

-.0010 


c2 

.9098 
.6452 

1.7743 
.1050 
.6322 

3.6292 
.17419 

-.0026 


"2 
.9661 
.6275 

1.7256 
.1765 
.6601 

3.8183 
.17286 

-.0049 


2t 


s  p 


If 

C, 


°2t  "  ° 
.8412 


It 


c2 

1.0194 
.6139 

1.6832 
.2416 
.6822 

3.9938 
.17081 

-.0083 


*2 
1.0700 

.6033 
1.6591 

.3016 

.6997 
4.1577 

.16828 
-.0116 


Example  IXb 

Solution  for  the  rotor  of  Fig.  5  (divergent  curved  channel 
walls)  with  a  radial  temperature  gradient  and  radial  coordinates  for 
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eq-aal  areas.     This  is  an  isentropic,  subsonic  solution  for  the  in- 
verse problem    (H(^)  given)  for  the  rotor.    One  complete  approxi- 
mation for  the  inverse  problem  for  the  rotor  is  given  below  to  il- 
lustrate the  iteration  procedure. 

Equations  (63),  (65),  (73),   (75),  (77)  (with  appropriate 
subscripts)  and  the  following  two  relations  are  used; 

>4 


s .  \(£lYJ>  Is  4&tLp)  MtiW.^l  '  (,M,kAdn  t  d.0, 


dffC,    (80) 


«/#-U*64 


(27). 


When  iterated  across  the  channel,  Equation  (80)  becomes: 


*  -tmif[^('-0^hH^^t]  -***> 


(81). 


+inmtmu^. 


Conditions  at  station  2  are  the  same  as  after  the  stator, 
Example  IXa  above.     The  work  output  function  Hfr£  )  is  given  in  this 
example,  and  the  blade  shape  (spouting  angle)  is  determined. 


Data  used  in  the  iteration: 

dLtzt 


m  =  . 029242 


<Hi-Hi-l>  =  0 


*& 


=  2.8593  H(/J  =  1.300<c-Snr 


given) 


cLH 
"6 


=  0 
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a. 


-6 


"(ft 
As 


d3     e3     £3 


y  ctGt&f-fJ,') 


fa  -<--,) 


7600   .8661    .9605    1.0464  1.1258   1.2000 


1.3000   1.3000   1.3000    1.3000   1.3000  1.3000 


1.5129   1.5626   1.6143    1.6671   1.7210   1.7752 


-.2129  -.2626   -.3143    -.3671   -.4210  -.4752 


-.2801   -.3032   -.3272    -.3508  -.3739   -.3960 


1.7889   1.7317    1.6766    1.6238  1.5736 


.86305  .86305    .86305    .86305  .86305 


-.02931    -.02801       -.02751      -.02681    -.02631 


2    ,* 


sffi=£) 


le     2 


.1061        .0944  .0859  .0794        .0742 


.01489      .01501         .01496        .01482     .01462 


2zLj-Uf\  .007820   .006658      .005821      .005185  .004682 
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Try  Cj  =.1039. 

a3       b3       c3       d3       e3       f3 
$3  .1039     .2379     .3295     .4159     .5242     .6891 

^3  -.2801    -.3032    -.3272    -.3508    -.3739    -.3960 

93  -.0156     .0188     .0932     .1948     .3354     .5513 

.1047     .1088     .1145     .1209     .1274 


A 


$3 

~93  .0016  .0560  .1440  .2651  .4433 


903 


.3242  .7881         1.1827        1.7707        2.9097 


"^T?  -.0005        -.0441         -.1703        -.4694      -1.2889 

^zr  -2.9845      -2.1545        -1.5785    -1.0580        -.4485        +.4290 

Tf(0)     +'4477        -1704  -.4466      -.4956        -.2870        +.3432 

^(f?*)      "-0260  '0409  -1804  .3698  .6316        1.0336 

^r  +.4217        -.1295        -.2662        -.1258        +.3446      +1.3768 


+.1461        -.1978        -.I960        +.1094        <-.  8607 


? 


903 


-(k^r  -.0065         +.0592         +.0700         -.0642         -.5647 

&Rr  f.  00076       +.00082       +.00008      -.00214      -.00809 
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The  constant  C?  is  found  by  applying  Equation  (77)  at  the  inner 
channel  boundary. 

^  =  3.  7894  -  . 14286     [. 0785  +  . C108  ♦  . 0002  +  2.  60o]  =  3. 4052 

=  .6878        C2  =(.6878  p~   =8986 


P3  =  .  5460    (1.0682)**-/ 
b 


7-' 


P3  r       .8986 


.6878 


'3 

.8994 

.6900 


"3 

.9002 

.6921 


a3 

.9003 

6924 


e3 

.8981 

.6865 


(-^■J      1.2597         1.1589         1.0947         1.0489         1.0063 


j. 


/-^M*  1.1793        1.1111         1.0668         1.0347         1.0045 
Av.  1.1452         1.0889         1.0507         1.0196         .9845 

V*3  *)'  1014  .1006  .0998 


.0997 


.  1019 


43 

.8900 

.6651 
•  9506 
.9645 

.  1100 


Av. 


.  1010 


.  1002 


.0998 


.  1008 


.  1059 


The  axial  velocity  integral  is  iterated  in  steps  across  the 
channel  as  follows: 

2. 0486  (. 08717  . 02931)  ♦  . 01550  =  . 13403 
1.  8856  (. 08648  -.  02801)  -  . 01867  =  . 09158 
1.7616  (.08613  -.02751)  -.01684  ■  .08642 
1.6556  (.  08700  -.02681)  +.00869  =  . 10834 
1. 5492  (. 09140  -  . 02631)  t. 06386  =  . 16470 

The  mass  flow  equation  is  iterated  to  solve  for  C,  as  follows; 

. 029242  =  . 017052  (. 06701  +  Cj)  t  . 016344  (. 17982  +  Cj) 
+  . 015718  (. 26882  t  Cj)  f  . 015110  (. 36620  +  Cj) 
*  .014393  (.50272  *  C^ 

cl       .0VS61?     " -1039 


<p3 

Average 

13403 

.06701 

22561 

.17982 

31203 

. 26882 

42037 

. 36620 

58507 

.  50272 
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Solution  at  station  3: 


i. 

.7600 

.8661 

.9605 

1.0464 

1. 1258 

1.2000 

<Ps 

.  1039 

.2379 

.3295 

.4159 

.5243 

.6890 

*3 

-.2301 

-.3032 

-.3272 

-.3508 

-. 3739 

-.3960 

63 

-. 0156 

.0188 

.0932 

.  1948 

.3354 

.5512 

P3 

.6878 

.6900 

.6921 

.6924 

.6865 

.6651 

r, 

3.4052 

3. 5688 

3. 7243 

3.8706 

4. 0010 

4.  0978 

°3 

.20197 

. 19335 

. 18583 

. 17S86 

.  17158 

. 16231 

G//,) 

4.619 

2.366 

1.922 

1.672 

1.434 

1.  167 

ft 

12°12' 

22°55' 

27°30! 

30°55' 

34°53» 

40°35' 

A  T3 

0 

+. 0064 

+.  0096 

+.0097 

+. 0070 

+  .0014 

Note:    A  supersonic  solution  is  indicated  for  problem  IXb  by  using 
average  values  in  the  mass  flow  and  energy  equations 
(  <f>    =    3.780).     However,   the  supersonic  solution  is  im- 
possible in  this  case  since  the  pressure  goes  to  zero  at  the 
outer  channel  boundary  before  the  solution  is  reached. 

Exai-n.pl  e  Xa. 

Examples  Xa,  b  are  the  isentropic,   subsonic  solutions  for  the 
stator  and  rotor  of  Fig.   5.     This  is  an  example  of  the  inverse  problem 
in  which  H{   £3     )    and        ^  =  0  are  specified  and  the  blade  shapes  of 
both  the  stator  and  rotor  are  thus  determined. 

Equations  (27),  (63  ),   (65),  (73  ),  (74),  (75)  (with  appropriate 
subscripts)  are  applicable  for  the  stator  in  addition  to  the  following 
relation: 
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■/, 


%--\\m^m^ 


d0. 


**  dfcAz) 


*Z      <Pe    **   *$■< 


d&i-C,         (82). 


Conditions  at  station  1  are  assumed  the  same  as  for  Example 


IXa. 


a2  °Z  c2  a2  e2  x2 

Mfs)  1.3000    1.3561    1.4061    1.4515    1.4935    1.5328 
(given) 


Ar      1-6561 


1.5962         1.5455         1.5024         1.4651 


1.4325 


m  =  . 029242 
dT,± 


*fi 


=  4. 0746 


«#- 


Solution  at  station  2: 


<P. 


%        .7850 

.8496 

.9098 

.9661 

1.0194 

1.0700 

4>&        .6791 

.6327 

.6004 

.5760 

.5532 

.5311 

/\r«i»0i.656l 

1. 5962 

1.  5455 

1.5024 

1.4651 

1.4325 

%       -• 0679 

.0228 

.0977 

.  U 

.  2177 

.2656 

P2        .  6367 

.6780 

.6932 

.7308 

.7489 

.7635 

Ta      2.3310 

3. 5501 

3.7485 

3.931C 

4. 1015 

4. 2627 

Q2        .19116 

. 19093 

. 18604 

.  18590 

.  18279 

. 17912 

F(Q  2.  439 

2.523 

2.574 

2.608 

2.648 

2.  697 

#»        22°18' 

21°36' 

21°15' 

20°59* 

20°41» 

20°21» 

&T2            o 

-.0035 

-.0052 

-.0051 

-.0038 

-. 0019 

Cj  • .6791         C2 

e  .8790 

P       -  P 
*2t         If 

Q2t  =  Clt, 

Z    -  T 
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Example  Xb. 

Solution  for  the  rotor  of  Fig.    5  for  the  inverse  problem  in 
which  H  (  £    )  and      h3=  0      are  specified  and  the  blade  shape   G  (fs) 
is  to  be  determined. 

Equations  (27),   (63),  (65),  (73),  (75),  (  77)  (with  appropriate 
subscripts  are  applicable  for  the  rotor  in  addition  to  the  following 


(Jk)  Ul  IjL  dZkt/hP7L  «L//  ^1  *  ** 


<*£+$ 


(83). 


Conditions  at  station  2  are  the  same  as  after  the  stator, 


Example  Xa. 


m  *  . 029242 

dttt         . 

8593          \3 

=  0     G<£ 

1   <h 

Solution  at  station  3: 

a3 

b, 
3 

c3 

d3 

e3 

f3 

£         > ?600 

.8661 

.9605 

1.0464 

1. 1258 

1.2000 

05        . 2985 

.3398 

.3627 

.3830 

.4110 

.4602 

9s      -.0048 

.0269 

.  1026 

.  1794 

.2630 

.3682 

)i3t          0 

0 

0 

0 

0 

0 

P3        . 6876 

.6873 

.6868 

.6857 

.6825 

.6758 

Q3          20196 

. 19280 

.18431 

.17763 

.  17105 

. 16416 

T3     3.4049 

3. 5639 

3.7388 

3.8602 

3.9940 

4. 1165 

G(jj)  2.  546 

2.549 

2.648 

2.732 

2.739 

2.608 

@3        21°25' 

21°23' 

20°41' 

20°07« 

20°02' 

20°57' 

&T3        o 

-.0001 

+.0011 

-.0010 

+.0001 

*.  0002 

Cl  = 

.2985 

C2  =  .8985 
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Example  XI 

This  example  has  been  constructed  to  illustrate  the  increase 
in  gross  thrust  obtainable  by  the  use  of  a  radial  total  temperature 
gradient. 

A  hypothetical  one-stage  turbine  consisting  of  a  stator  and  a 
rotor  has  been  assumed  similar  to  that  shown  in  Fig.   3.     In  order 
to  give  this  example  more  physical  meaning,   rotor  blades  have  been 
used  which  are  tapered  linearly  in  area  in  such  a  way  that  the  ratio 
of  hub  to  tip  cross -sectional  area  is  three.     Using  a  rotor  tip  speed 
of  1100  feet  per  second,  the  stress  due  to  centrifugal  forces  was 
calculated  at  various  points  along  the  blade.     Using  the  high  -tempera- 
ture data  for  "Vitallium",  (H.£.    21)  (NR-1G),   given  in  References  16 
and  17,  it  was  then  estimated  that  the  blade  root  could  safely  stand  a 
temperature  of  1450  F,  while  the  blade  tip  temperature  could  be 
1850°F. 

Using  this  data,  the  thrust  of  this  turbine  is  calculated  for 
two  conditions:    (a)    A  constant  radial  total  temperature  of  1450  F, 
and  (b)  a  linear  temperature  gradient  of  400  F  with  the  same  root 
total  temperature  of  1450  F.     The  mass  flow,  boundary  configuration, 
combustion  chamber  total  pressure,  and  work  extracted  arc  the  same 
in  both  cases. 

Exam.pl  e  XIa. 

Solution  for  the  stator  and  rotor  of  Pig.    3  with  no  radial  total 
temperature  gradient  and  equal  area  coordinates.     This  is  an  isen- 
tropic,   subsonic  solution  to  the  inverse  problem  (H    (  f3)      given). 
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For  no  radial  total  temperature  gradient,  H  (    $      )  =  /  J 
=  constant,  and  a  cylindrical  channel,  Equation  (82)  reduces  to 
4>z  -  C.  .     Equations  (63),  (64),  (65)  and  (67)  are  also  U6ed  in 
the  solution  for  the  stator  . 
Conditions  at  station  1:    (given) 


£*  =  2.7071 
2/  =  2.6726 
u    =  880  ri 

cj/P  =  1100 


Flt  =  l.00 


Cu  =  .36940 


Pj     =     .9562      Qj    =   .35775 


<£,=  .4915 

A/S  o 


sec 

ft 


Tu  =    1910°R  (const.)  W  =  148.26 


sec 


sec 
/?„  =  1.25  ft 


p..    =    16,300  -1^- 
It  -£T 

m       =  .031650 


H(£)=  1.3964  (const.) 


\3=  0 


Solution  at  station  2: 


i 


2 


"2 
.8000 

.7260 
1.7455 

.4812 
C2  .21909 

T2  2.1965 

F  =-~   2.404 


22°34' 


2 


1.071 
0 


u2 
.8544 

.7260 

1  .6342 

.5236 

.23268 

2.2501 

2.251 


23°59' 


1  .005 
+.0002 


fc2 

.9055 

.7260 

1 . 542 1 

.5585 

.24369 

2.2921 

2.124 


25°13' 


.951 
0 


u2 
.9539 

.7260 

1.4639 

.5879 

.25277 

2.3259 

2.016 


26°25» 


.904 
-.0003 


c2 

1.0000 

.7260 

1.3964 

.6125 

.26027 

2.3533 

1.923 


27°30' 


.873 

-.0001 


Cj  =  .7620 


?*  =  ?/t    j     P2t  =  *  It'  °2t        "it 


C2  =  .811- 
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For  the  rotor  with  no  radial  total  temperature  gradient, 
H      if*)  s  constant,   and  A^     =0,   the  solution  at  station  3 

is  easily  found  to  be: 


*i 


*  3 
2i 


■       .7806 

0 
=      .5003 
=    2.2211 


C3t 
?3t 


=      .22525 
=  2.  3081 
=      . 24798 
=      . 5724 


G  = 


fr 


.    4, 


1.025        1.095 


44°20' 


42°25' 


1.160  1.222       1.281 

40°46'        39°17'      37°56' 


Assuming  isentropic  expansion  to  atmospheric  pressure  through 
an  ideal  nozzle,  the  total  gross  thrust  is  found  to  be  11 ,  966  pounds. 
Example  Xlb. 

ution  for  the  same  stator  and  rotor  as  Example  XIa  except 
with  a  linear  radial  total  temperature  gradient  of  400  F .     The  mass 
flow,   combustion  chamber  total  pressure,  and  work  extracted  are  the 
same  as  in  the  above  example. 

Equations    (27),  (63),  (64),  (65),  (67),       and  ( 82 )  are  applicable. 

Conditions  at  station  1:    (given) 


t 


Q 


It 


°i 


.8000 
2.7071 
1.00 
. 36940 
.9512 
2. 6687 
. 35643 
.5185 
0 


"1  wl 

.8544        .9055 
2.8614     3.0062 


1.00 
. 34948 
.9512 

2.8208 
.33721 
.5331 
0 


1.00 
.33265 
.9512 

2.9635 
. 32097 
.5467 
0 


"1 

•  9539 
3. 1434 
1.  00 
.31813 
.9512 
30988 
.30696 
.  5587 
0 


cl 

1.0000 

3.2741 

1.00 
. 30543 
.9512 

3.2276 
.29471 
.5705 
0 
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m  =  .031650       ^r*~=2.  8350   H  (4 )     =  1.  3964  (const  .),  A  ■ 

CLff  3  3 


Solution  at  station  2 : 


a2 

b2 

C2 

d2 

e2 

*. 

.5844 

.6975 

.7955 

.8819 

.9589 

A2 

1.7455 

1.6342 

1.5421 

1.4639 

1.3964 

P2 

.  5018 

.5443 

.  5774 

.6036 

.6249 

U 

2.2231 

2.4050 

2.5697 

2.7212 

2.8625 

°2 

. 22574 

.22632 

.  22474 

.22183 

.21829 

^2 

2.987 

2.342 

1.938 

1.659 

1 .  467 

ft 

18°30' 

23°07' 

27°  18' 

31°09' 

34°19' 

M2 

1.044 

.  974 

.920 

.878 

.850 

*£ 

-. 0001 

+.  0054 

h  0064 

+.  0049 

+  .0017 

Cl  = 

.  5844 

r*  = 

Ttt    '  P2t 

s  Plt'   Q2t 

it 

C2^ 

.8212 

For  the  rotor  with  H  (£3)  =  constant  and       A^    -  0, 

Equation  (64)  reduces  to  P-  «  constant.     Equations  (63),  (65),  (71), 
and  the  following  integral  are  applicable: 


■k 


Solution  at  station  3: 
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<t>3 

°3 
P 


3t 


3t 

4>3 


.6002 
0 

.5289 
2.2566 

.23435 
2.3081 

.5724 

.24798 
1.  333 


36°52 


.7426 
0 

.5289 
2.3853 

.22171 
2.4624 

.5912 

.24005 
1.151 


40°58 


-.0017 


.8514 
0 

.5289 
2.5060 

.21103 
2.6072 

.6075 

.23300 
1.064 


43°13' 


-.0023 


.9397 
0 

.5289 
2.6203 

.20182 
2.7444 

.6220 

.22658 
1.015 


44°35' 


-.0020 


1.0141 
0 
.5289 

2.7293 
.19376 

2.8751 
.6345 
.22067 
.986 
45°22' 

-.0011 


Cx  =   .6002 


C2  =  .8336 


Assuming  isentropic  expansion  to  atmospheric  pressure 
through  an  ideal  nozzle,  and  summing  the  thrust  by  stream  tubes 
the  total  gross  thrust  is  found  to  be  12,883.8  pounds.    This  is 
an  increase  of  917.8  pounds,  or  7.7  percent,  over  than  found  in 
Example  XIa. 
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Mi.  *• 

a 

.8000 

b 

.8544 

c 

.9055 

d 

.9539 

e 

1 .0000 

e,/  //////  /   *>/  /   /  /  /  /  /e,/    // 


a> 

o 

c 


o 
o 
o 


"g    0.70 


0.60 


O.IO  0.20  0.30 

$ 
Axial     Coordinate 


0.40 


0.50 


FIG.3—  CHANNEL    CONFIGURATION     AND    COORDINATES     FOR 
EXAMPLES     I,  H,  AND  31 
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*, 

h 

*3 

/Jf) 

a 

.800 

.790 

.780 

-.050 

b 

.850 

.850 

.850 

0 

c 

.900 

.910 

.920 

+  .050 

d 

.950 

.970 

.990 

+  .100 

e 

1  .000 

1  .030 

1  .060 

+  .150 

1.10 


O.IO 


0.20 


0.30 


0.40 


0.50 


Axial    Coordinate 


FIG. 4  —  CHANNEL     CONFIGURATION    AND    COORDINATES    FOR 
EXAMPLE    HI 
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1.30 


0.10 


0.30 


0.20 
Axial    Coordinate 


0.40 


0.50 


FIG. 5  -  CHANNEL    CONFIGURATION    AND    COORDINATES    FOR 
EXAMPLES  ZEE  -31 
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Examples   JZ  ,  X 


5, 

Si 

iz 

^3 

(JL) 

(81,3 

a 

.800 

.785 

.760 

-.05 

-.10 

-.15 

-.25 

b 

.840 

.842 

.848 

0 

.02 

.04 

.10 

c 

.880 

.899 

.936 

.05 

.14 

.23 

.45 

d 

.920 

.956 

1  .024 

.10 

.26 

.42 

.80 

e 

.960 

1.013 

1  .112 

.15 

.38 

.61 

1  .15 

f 

1  .000 

1.070 

1  .200 

.20 

.50 

.80 

1  .  51 

Examples  21  -X 


V, 

*. 

43 

VT^l 

(!U 

\d&)3 

U  i% 

(&) 

u% 

a 

.8000 

.7850 

.7600 

- .0500 

-.1000 

-.1500 

-.2500 

-.2500 

-.2500 

b 

.8438 

.8496 

.8661 

.0047 

.0360 

.0791 

.1332 

.1467 

.1720 

c 

.8854 

.9(98 

.9605 

,0567 

U1627 

.2829 

.4972 

.5163 

.5474 

d 

.9252 

.9661 

1 .0464 

.1065 

.2813 

.4684 

.8455 

.8620 

.8891 

e 

.9633 

1 .0194 

1  .12  58 

.1541 

.3935 

.6398 

1 .1789 

1.1893 

49 

f 

1 .0000 

1 .0700 

1  .2000 

.2000 

.5000 

.8000 

1.500C 

1.5C00 

1 .5000 

TABLE  X  -  STREAMSURFACE    COORDINATES,  SLOPES   AND   CURVATURES 
USED    FOR    EXAMPLES  IZ -X  AND  FIGURE  5 


5 
o 


to 
to 


Corrected  |  2 

Prob.IZ  b 

Prob.lL  b 

© 

a 

.785 

.785 

0.028 

b 

.  82  0 

.838 

c 

.866 

.90  3 

d 

.923 

.968 

e 

.990 

1  .026 

O  OP4 

/ 

\y.v/C~ 

f 

1  .070 

1  .070 

0.020 

/     / 

m(l2) 

Station  27Prob.32" 

a/ 

/     i         j 

0.016 

1      (Su 

Dsonu 

:; 

/        A    -"(^ 

'              /Station3,Prob. 

3Tb 
:) 

/  A 

(Supersonic 

Fl 

/                  \ 

0.012 

o 

Sta 

Vm(    §3) 

tion3,Prob.:DZb~ 
(Subsonic) 

0.008 

Ay 

/  / 

/ 

r\  r\dA. 

U.UUf 

n     r 

1 

0.76 


0.84 


0.92 


4 


.00 


1.08 


16 


Radial    Coordinate 

FIG.6-MASS    FLOW    FUNCTION     FOR    SUBSONIC    AND    SUPERSONIC 
ROTOR    FROM    SUBSONIC    STATOR 


1.20 


99 

Streomiines  Originally  Assumed 
Corrected  Streamlines, Prob. 3Tb 
(Supersonic  Mg) 

Corrected  Streamlines, Prob.JZ'b 
(Subsonic  Mg) 


0.80 


FIG.  7  -STREAMLINE     SHIFT  THROUGH    ROTOR 
FOR   SUBSONIC  AND    SUPERSONIC 
"GOVERNING"  MACH   NUMBER 


0.75 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


Axial    Coordinate 
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2  2   ^S 


V2  +  A2 

3  3 


Stator  :  Cot  A  (  «2)  =  F(  l2)  =  -^-        Rotor:  Got  /*,(!,)  =  G(  L)~-  ^3 


2'       0 


3W3' 


<*: 


FIG.  8   -TYPICAL    TURBINE    VELOCITY    DIAGRAM     SHOWING    SIGN 
CONVENTION    USED    IN    ALL    EXAMPLES 


y2  2 


V^2>A2 
r3  3 


Stator:  Cot/?2(  §  2)  =  F(  4  2)  =  -^r~        RotorCot/53(  ^3)=6(l3)=^— 

9  2 


FIG. 9    -TYPICAL     COMPRESSOR     VELOCITY    DIAGRAM    SHOWING 
SIGN    CONVENTION 
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